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Abstract 

This work reports on the experimental investigation of surface plasmon polaritons (SPP) on 
cylindrical wires of small diameters. It has been accomplished within the scope of a diploma 
thesis at the Max-Planck Research Group for Optics, Information and Photonics of the Friedrich- 
Alexandcr University of Erlangen-Nuremberg. 

Applying a new technique that was developed by this group and recently reported 7^, 84, 



l63l |82| , single wire and wire array samples of gold (Au) and silver (Ag) with small diameters 
(400 nm < < 3 urn) and high aspect ratios (< 75 000) in photonic crystal fibers and single 
hole capillaries arc fabricated. Additionally, effective bulk metal wires of Au and Ag are created 
by a hybrid technique, including fiber tapering and magnetron sputter deposition for a large 
number of different diameters between 13 urn to 50_um. 

First results of optical in-fibcr excitation of SPP [84J suggested an in depth study of the disper- 
sion relation of surface plasmon polaritons. The measurement of the absorption, corresponding 
to the excitation of SPP modes was realized with a new goniometer-based experimental setup 
for oblique incidence laser beam scattering at a single wavelength. A device was developed that 
allows a simultaneous measurement of the scattering amplitude in s- and p-polarization under 
varied incident angle. 

For comparison, a model [7(| for the dispersion of discrete SPP modes and their angles of 
excitation is derived. The computed results are for the first time compared with a theoretical 
approach to electromagnetic scattering theory for cylinders, that is calculated analogue to Mie 
scattering [H^ and to an exact solution of Maxwell's equations for cylindrical SPP. The validity 
of the model is proved for higher incident angles. The first oblique incidence examination of 
SPP on metal wires for this size range is reported. The experimental results for low incident 
angles show good accordance to the numerical results from scattering theory. 



Zusammenfassung 



Das Thema dieser Arbeit ist die expcrimcntcllc Untcrsuchung von Plasmon-Polaritonen auf dcr 
Obcrflachc von Drahten geringen Durchmcsscrs (SPP). Die Untcrsuchungcn wurden im Rahman 
einer Diplomarbeit an dcr Max-Planck-Forschungsgruppe fur Optik, Information und Photonik 
der Friedrich-Alexander-Universtat Erlangen-Niirnberg durchgefiihrt. 

Auf Grundlage einer neuen Technik, die an dieser Gruppe entwickelt wurde (7^, 84, 63, 82j, 



werden Einzcldrahtc und Drahtarrays aus Gold (Au) und Silber (Ag) mit geringen Durchmessern 
(400 nm < < 3 |J.m) und grofiem Langenverhaltnis (< 75 000) mittels Matrizen gezielt gefer- 
tigter Photonischer Kristallfasern hergestellt. Aufierdem wird eine kombiniertc Technik aus 
Fasertapern und Magnctronsputtcrn zur Hcrstellung von Probcn mit Durchmessern zwischen 
13 urn und 50 |am angewandt. 

Erste Ergebnisse aus der Anregung von SPP durch glasfasergefiihrtes Licht [84| legten eine de- 
taillicrterc Untcrsuchung dcr Dispcrsionsrclationcn solchcr angeregten Obcrflachcnmodcn nahc. 
Die hier angewandtc Technik, diese anzuregen, basiert auf einem neuen, goniomctcrbasicrtcn 
Aufbau zur Messung der Streuamplitudc cincs schrag zur Drahtachsc cinfallcndcn Laserstrahls 
bei konstanter Wcllcnlangc. Das optische System wurde so entwickelt, dass eine simultanc 
Messung der Streuamplitudcn cinfallcnder s- und p-Polarisation unter Winkclvariation moglich 
ist. 

Ergebnisse auf Basis eines Modells [84| soldier SPP-Moden werden erstmals mit numeri- 
schen Ergcbnissen verglichen, deren Berechnung nach einer Hcrlcitung von winkclabhangigen 
Streuamplituden und Feldverteilungen fur die Streuung an cincm Zylindcr analog der Mie- 
Losung [56j der Maxwellgleichungen erfolgt. Aufierdem erfolgt ein Vergleich mit der exakten 
Losung dcr Maxwell- Glcichungcn fiir SPP. Dcr Giiltigkcitsbcrcich des Modells wird dabci auf 
groficrc Einfallswinkcl bestimmt. Hicrmit wird erstmals die Anregung von SPP auf Drahten 
dieser geringen Grofic untcr schragcr Inzidenz untersucht. Die Messergebnissc fiir geringen 
Einfallswinkcl zeigen gute Ubercinstimmung mit den Bercchnungcn der Streu-Theorie. 



Contents 



9 



Contents 



1. Introduction [l 

1.1. The examination of surface plasmon modes on wires Q 

1.2. Overview and structure of the thesis [i 

2. Theoretical background jj 

2.1. From Maxwell's equations to surface plasmons 5_ 

2.1.1. Maxwell's equations 5_ 

2.1.2. Light-matter interaction g 

2.1.3. Dispersion |8 

2.1.4. The vector and scalar wave equations ]J 

2.2. Planar surface plasmon polaritons 1J 

2.2.1. Background 1J 

2.2.2. Dispersion of planar surface plasmon polaritons 2£ 

2.3. Surface plasmon polaritons on a wire 25. 

2.3.1. Background and validity 25. 

2.3.2. Model for the dispersion relation 2J 

2.4. Scattering of light by a small cylinder 2J 

2.4.1. Background 2J 

2.4.2. Cylindrical coordinate system 21 

2.4.3. Field expansion in vector cylindrical harmonics 31 

2.4.4. Boundary conditions and analytical solution 3J 

2.4.5. Perpendicular incidence 3J 

2.4.6. Scattering efficiencies 36 

2.5. Asymptotic scattered field and intensity pattern |41 

3. Sample fabrication 

3.1. Material properties 



3.2. Free-standing wires 

3.2.1. Tapering 

3.2.2. Sputtering 

3.3. Wires in fiber silica matrix 

3.3.1. Fiber drawing 

3.3.2. Fabrication of wires in fibers 

4. Experimental setup 

4.1. General setup 

4.2. Beam optics 

4.2.1. The light source 

4.2.2. Polarization labelling by beam chopping 

4.3. Scattering 

4.3.1. Optics 

4.3.2. Free-standing scattering setup 

4.3.3. In-liquid setup 

4.3.4. Visual process observation 

4.4. Measurement and referencing 

4.4.1. Optical detection 

4.4.2. Principles of lock-in amplification 

4.4.3. Polarization dependent lock-in detection and noise handling 

4.5. Data acquisition and control software 

5. Theoretical results 

5.1. Numerical computation of the dispersion model 

5.2. Comparison to an exact analytic solution 

5.3. Numerical computation of the scattering coefficients 

5.3.1. Inherence of SPP resonances 

5.4. Numerical computation of the scattered field distribution 

6. Experimental results and discussion 

6.1. Observed scattering pattern 

6.2. Measurement of resonance angles 

7. Conclusion and outlook 

A. Mathematical and electromagnetic foundations 

A.l. The Bessel differential equation and its solutions 

A. 2. The Poynting vector of plane waves 



1 



1 



Chapter 

Introduction 

1.1. The examination of surface plasmon modes on wires 

Surface plasmons, collective vibrations of the electron gas propagating on the surface of 
metals, have experienced a renaissance during the last years since their first experimen- 
tal exam inati on and theoretical treatment in modern physics by Sommerfeld [92] (1899), 



Zenneck [112|] (1907) and Mie [5Q] (1908). 

The optical effects of surface plasmon (SP) excitations, particularly the colour effects 
of metal colloids were phenomenologically known by craftsmen much earlier. One of the 
first and most impressive applications of dichroic glass is the 'Lycurgus cup' from the late 
roman era (4th century AD), which shows a different colour for transmitted and reflected 
white light. 

Today, since the initial experimental studies on Raether [73| , SP are finding an increasing 
number of applications 0. Especially the advance in optical technology recently allows 
an in depth optical examination of surface plasmons that were before in a wide range of 
frequencies preferentially examined by electron diffraction and loss measurements. 

The advantages of surface plasmons can be attributed to its intrinsic nature of being 
confined in the near vicinity of metal-dielectric interfaces. This confinement leads to a 
strong enhancement of the electromagnetic fields, which makes SP extraordinarily sensitive 



to surface conditions. This sensitivity can be utilized e.g. for sensing applications [65| 



Surface plasmons also offer new possibilities for nanoscale photonic circuits [a], for the 



fabrication of subwavelength waveguides [ljj and may offer possibilities to realize a basis 



for optical switches and circuits below the diffraction limit of the applied wavelength [41 1 . 

In addition they are shown to be well applicable to construct collectors for electromag- 
netic radiation near the surface of samples. These allow the measurement of near field 
distributions and reveal new insights into the radiative properties of samples, scanning 
near field optical microscopy (SNOM), an application that has already entered commercial 



2 



1. Introduction 



development. 

The progress lead finally to the concept of specifically 'designed materials' to promote 
desired dispersions that are not achieved with natural materials [661 ] . 

The common focus of the theoretical examination of surface plasmon polaritons (SPPs) 
always was on spherical particles as these are of special importance in many interdisci- 
plinary applications from biology, the treatment of organic cells to atmospheric physics 
and aerosol science. 

Nevertheless, the theoretical examination of surface plasmon modes on cylinders over 
the last few years gathered increased attention throughout literature [ 



49|, [85|, 123, 181| 



llll . l62l . l85| |. Experimental studies for many of these calculations still have to be ac- 



complished, although in some geometries first results are promising for further in depth 
examination [25l, 111 . 



Recently by Schmidt et al. 84J the fabrication and examination of metal wire arrays with 
an unsurpassed high length to diameter ratio and good surface smoothness was reported. 
These samples were fabricated by a new technique inside the hole structure of photonic 
crystal fibers (PCF) @. 

The resonance frequencies of SPPs on these wires were determined by measuring trans- 
mission spectra of the fibers. The results show good accordance with theoretical multipole 
calculations and a proposed 58, 84, 83 1 model for the description of the dispersion of 
SPPs. 



The aim of this work is to answer a number of questions that appeared in the course of 
these studies and are not yet treated by any other experimental work: 

1 . What is the range of validity of the proposed model for the propagation of radiative 
SPPs on metal wires? 

2. Does a description of the scattering of light on metal wires comprise SPP resonances? 

3. Does an experimental scattering examination of exciting SPPs on small metal wires 
reveal resonances that match the predicted dispersion relation in a range of wave 
vectors that is not feasible by in-fiber coupling? 

4. Is it possible to relate the observations from scattering to a near field description? 



1.2. Overview and structure of the thesis 

This thesis comprises the following parts: 

As a background for the further theory, in the first part of the thesis the theoretical de- 
scription of electromagnetic waves is introduced in general and the linear optical properties 
of matter, especially for metals are discussed (sec. [2]). 
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Subsequently (sec. l2.2[ ) an overview over the occurrence of planar surface plasmon po- 
laritons is given. This theory is then (sec. I2.3."2"j) applied to the special case of surface 
plasmon polaritons on a closed two dimensional surface, leading to the dispersion relations 
for discrete surface plasmon modes on cylindrically curved surfaces. 

From electromagnetic scattering theory (sec. I2.4|) an analytic approach, analogue to the 
Mie solution j5^j for infinite right circular cylinders is derived. 

The range of validity of the model is evaluated by a comparison of computed results 
for the applied materials and experimental parameters (sec. 15. l|) to the exact solution of 
Maxwell's equations for radiative SPP modes in the same case (sec. 15. 2| ). 

From the derived Mie theory, scattering field distributions and absorption patterns are 
computed and compared to the exact solution and the model (sec. I5.3|) . An analysis of the 
near field distributions and an approach for their explanation and interconnection to the 
mode excitation is presented. 

To test the theory, two different types of samples are fabricated, 

1. freestanding, metal coated silica wires for diameters D > A (sec. I3.2p and 

2. nanowires embedded in silica fiber for diameters D < A (sec. I3.3[) . 

The functional principle of the experimental measurement setup for probing these sam- 
ples is explained (sec. 0J). Results of the accomplished experiments are presented and 
discussed in section [6] with a comparison to the derived theories. The conclusion from the 
experimental and theoretical investigations closes the work (sec. [7} and propositions for 
further studies of the observed phenomena are given. 
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Theoretical background 



"Whether you can observe a thing or not depends 
on the theory which you use. It is the theory 
which decides what can be observed." 
(Albert Einstein, 1926) 

2.1. From Maxwell's equations to surface plasmons 
2.1.1. Maxwell's equations 

Electromagnetic waves are in general described by equations, which were first collated from 
previously already known physical laws to a consistent system by James Clerk Maxwell in 
1865 jEHj. Usually they are written in the modern form of four lineai0 differential equations 
B 

X7-D = p (2.1a) 
d 

VxH-—D=j (2.1b) 
d 

VxE + —B = (2.1c) 
dt 

X7-B = (2. Id) 



A list of mathematical and physical symbols and abbreviations that are used in the following chapters 
can be found in appendix IB . 1 1 with references to the page where they are first introduced. 

*As the equations are linear, linear combinations of their solutions solve the equations again. 

2 In theoretical physics often the cgs system is applied. Here for a better insight into the experimental 
meaning of the results all derivations base on the SI notation. For an in depth discussion it is referred 
to \M 
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Coulomb's law (I2.1a|) expresses that the source for the electric displacement vector D is 
the free electric charge density p. The extended Ampere's law (12.1bp states that the vortices 
of the magnetic vector field H as well as change in an electric vector field is caused by the 
electric current density vector j. The vortices of the electric vector field E are produced by 
a change in the magnetic induction vector B and vice versa, following Faraday's law fj2. lc() . 
Gaufi' law states that the magnetic field is solenoidal; there are no magnetic monopoles 
(|2~Tdf l. 

In vacuum the relations of electric displacement and electric field to magnetic induc- 
tion and magnetic field are connected by the dielectric permittivity of the vacuum £q = 
8.8542 x 10~ 12 As/ (Vm) and the magnetic permeability of the vacuum fio = 47rxl(r 7 Vs/(Am): 



D = e E B = hqH (2.2) 

The vacuum speed of light cq = 2.997924 58 x 10 s m/s in the SI system defines a basic 
physical constant which can be determined by measurement. From a measurement of £o> 
Ho can be calculated using the relation: 

1 

Co = ; 



2.1.2. Light-matter interaction 

It is possible to introduce the response of matter to electromagnetic radiation by introduc- 
ing additional terms into equations 12.21 called the electric polarization P and the magne- 
tization M: 



D = e E + P B = po(H + M) 

Since the wavelength of visible light is large compared to atomic distances, it is often 
sufficient to treat the response of media to light macroscopically. The electric polariza- 
tion P and the magnetization M are defined as functions of the electric field P(E) and 
the magnetic field M(H), representing the material response. This dependency can be 
expressed in good approximation, applying an expansion in power series in terms of sums 

3 No magnetic monopoles have been verified up to now. Indeed, an elegant theoretical explanation of the 
origin of the quantization of electric charges by Dirac JTi|] is based on the assumption that magnetic 
monopoles actually exist. An idea that is well described in [3l| 
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over vector entries: 

3=1 j=l k=l j=l k=l 1=1 

= Mo + 1 + EE + xgLw*. + ■ • • 

j=i j=i fe=i i=i fc=i i=i 

Thus the dependency is written as a successive sum over dielectric susceptibility tensors 
T]j2 1 and magnetic susceptibility tensors Xjl i P^ us a ^° ias OI " polarization Pq and mag- 
netization Mq. The terms of n = 1 are responsible for a linear dependency on the electric 
and magnetic field. Higher order terms (n > 1) describe nonlinear effects, which can be 
neglected as well as the bias terms for most optical materials under the assumption of field 
strengths which are reasonably lower than inner atomic Coulomb fields. 

Due to the isotropic nature of the materials that are treated throughout the course of 
this work@, the two dimensional 3x3 tensors r] and X can be written as scalar quantities. 



P = rjE (2.4) 
M = X H (2.5) 

Assuming a homogeneous material, both are additionally independent of the spatial 
coordinates. With those material constants the dependency of D(E) and B(H) for elec- 
tromagnetic waves in linear, homogeneous and isotropic media can be defineqj- 

From equations (|2.41 12. 5[) the material dependent dielectric permittivity e and the mag- 
netic permeability tensor respectively scalar \i are defined: 



e:=7]/e + l \i := % + 1 

D = e$eE B = fj,Qfj,H 

Different to the dielectric permittivity e > 1, the magnetic permeability in physically 
real cases is only limited to be positive < fx as the magnetic susceptibility complies 
— 1 < x p 0- This is the mathematical expression of para- and diamagnetism respectively. 
Only by so called synthetic meta-materials, this rule is broken, which also allow classically 
forbidden negative refractive indices. Quantitatively for all non ferromagnetic substances 

4 Otherwise optical effects like birefringence result (for further details see e.g. [iB], ch. xi] or [l^]). 

5 It shall be pointed out that throughout this work e and /x are defined as relative, dimensionless values. 

For details it is referred for appendix lB.il 
6 See 0, §31] for details. 



8 



2. Theoretical background 



the magnetic permeability can be assumed to be unity in most cases in optics as \[i — 1 1 < 
10~ 6 . 

The material properties also influence the current density which is generally not zero 
but j = crE. This relation is also known as the differential form of Ohm's law, where 
a represents the specific conductivity of a material. The value of this material constant 
distinguishes conductors a ^ 0, among which metals show extraordinary high a, from 
insulators a 0. 

By the so called 'Maxwell relation' the complex refractive index can be defined from e 
and (i as an alternative quantity for the material response to electromagnetic waves: 



Its real part n r = cn/c represents the quotient of the vacuum speed of light over the 
phase velocity in the specific material. Its imaginary part k, which is also called 'extinction 
coefficient' represents the absorption in the medium, as described by Beer's law (see p. H7|) . 

2.1.3. Dispersion 

For the introduction of the dielectric susceptibility the material properties were described 
macroscopically. This approach leads to correct results as it equals a spatial summation 
over the response to the incident field of a large number (N > 1 x 10 6 jl| of electrons and 
atomic nuclei. This number is sufficient to smooth microscopic fluctuations which allows 
a macroscopic treatment in the nonretarded limit for static fields or field fluctuations that 
are slow enough for the material to respond instantaneously. 

Indeed, propagating electromagnetic fields intrinsically cannot be described electrostati- 
cally. All media aside from vacuum show a dependency of the optical material properties 
on the frequency of the electromagnetic field. This dependency is called dispersion. Nev- 
ertheless, with this effect under consideration all derived relations remain valid, only the 
dielectric permittivity and magnetic permeability have to be described as frequency depen- 
dent quantities e(u>), 

The reason for dispersion is obvious from the point of view of a microscopic description 
of the material. For an e.g. harmonic time dependency an additional factor enters the 
electromagnetic fields E = E$e~ lult . Including this temporal dependency, the material re- 
sponse can no longer be assumed to be instantaneous, a phase discrepancy between driving 

7 An approximate number of the atoms inside a macroscopical volume with the size of optical wavelengths. 

A concise microscopic derivation can be found in [3lL Ch. 6.6]. 
8 To be precise, this approach is only valid as long as monochromatic waves are treated, otherwise it is an 

approximation. 




(2.6) 
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field and radiated field enters. For a quantitative derivation a microscopic description of 
the material is necessar^. 

The polarization P (eq. \2A\ of a volume element of a medium, can therefore be written 
as a superposition of molecular polarizations, expressed as the spatial mean molecular 
dipole moment (p) of N molecules with a molecular polarizability 7 mo i in that volume: 



P = N(p) 

{p) = £o7mol(# + Ei) 

Generally not only the macroscopic field E causes the molecular polarization p but 
intermolecular interactions add a second term E\. Applying the local field approximation^!, 
a relation between the polarization and the dielectric susceptibility is obtained: 



P = N 7mol (e E + 1/3P) 
El Nj mol 



n 



1 - 1/3JV7, 



mol 



The molecular polarizability 7 mo i in this expression represents the same information 
about the material response as the dielectric permittivity. The equation that connects 
7 mo i and e is called Clausius-Mosotti relational and includes the density of the material by 
the number iV of molecules inside the considered volume: 

3 (e- 1 

7mol 



N \e + 2 / 

For media of low density however the neighbour interaction is negligible {E\ ~ 0). 

A simple oscillator model is a sufficient approach for many applications to derive the 
frequency dependency of the dielectric permittivity. The electrons of all molecules are then 
assumed to be bound, each in a Lorentz oscillator potential. This adiabatic approximation 
bases on the assumption that the motion of the nuclei can be neglected. In fact, this is 
true up to very long wavelengths (infrared) due to the very different masses of electrons 
and atomic cores- 12 1. 



9 The dispersion is derived classically, omitting quantum mechanical effects, as a thoroughly quantum 
mechanical approach follows the same derivation with a replacement of the classical values by quantum 
mechanical expectation values and leads to the same results. For a derivation that accounts for quantum 
mechanics it is referred to Becker || and Adler 
10 See [U, P- 186] for details. 

11 For the assumption of optical frequencies, therefore fi « 1 and n 2 = e/eo this equation is also called 

Lorenz-Lorentz relation. 

12 Detailed treatment of the Lorentz oscillator and the free electron gas approximations can be found in the 

classical solid state literature, e.g. 0, 
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Figure 2.1.: Real Re(e) and imaginary part Im(e) of the dielectric permittivity at a resonance 
frequency ojq in the Lorentz oscillator model, including damping 7 7^ oF^I 

The equation of motion of an electron with the charge — e, mass m e and displacement 
r from its equilibrium position is an inhomogeneous linear differential equation of second 
order: 



—eE(r, t) = m c (r + jr + oj^r) 

with the resonance frequency ojq, considering a damping coefficient 7 

Solving this equation for the assumed time harmonic field E = Eoe~ lujt results in a 
molecular dipole moment p and via a summation to the macroscopic polarization P of: 

e 2 E 

P = Np = N 



m e uJq — lo 2 — itoj 



Applying this equation to equation (|2.4p the explicit dependency of the dielectric per- 
mittivity s(uj) = ?7/£o + 1 on frequency is obtained: 

N e 2 1 

e{u) = l + - — - ^— (2.7) 



13 The dispersion curve is calculated with numeric values that are not representative for real materials but 
that point out the resonance behaviour. A plot of the real dispersion for silica glass from the Sellmeier 



equation is depicted in figure 1371 
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Neglecting damping (7 = 0) the dielectric permittivity becomes infinite at the resonance 
frequency (u = ujq). In fact, damping cannot be omitted as it accounts for all dissipative 
processes (e.g. due to collisions between molecules) occur. Including these damping effects 
(7 7^ 0) the polarization becomes complex. 

The imaginary part lra(e(uj)), as displayed in fig. 12.11 is strongly peaked at a value of uj 
that is slightly smaller than the resonance frequency, indicating maximum excitation, while 
the characteristic of the real part Re(e(u;)) is slightly more complicated. For low frequency 
uj <C ujo the DC property of the medium prevails. 

Approaching the resonance frequency uj < ujq the Re(e(u;)) increases and reaches a 
maximum at uj = ujq — 7/2. This positive gradient behaviour is called normal dispersion 
(ND) and is the characteristic of most transparent substances in the visible wavelength 
regime. In this regime e re d < £biuej therefore different spectral components exhibit different 
phase velocities c{uj) and separate as soon as they pass material interfaces that are not 
parallel. For further increasing frequency ujq — 7/2 < uj a negative gradient follows, a 
region of so called anomalous dispersion until the frequency approaches uj = ujq + 7/2 and 
a second normal dispersion region begins. 

For frequencies uj ujq (ac, i.e. high frequency limit) damping becomes negligible, e{uj) 
becomes predominantly real and equation (12. T[) can be significantly simplified to: 



« 1 - 4. "I = — ' ( 2 - 8 ) 

uj 2 p e m e 

where oj p represents the plasma frequency of the material. 

For real materials the previous assumption of a Lorentz resonator at resonance frequency 
ujq must be extended. It is a good approach to consider a model in which iV molecules, each 
with Z electrons are situated in a volume element. Of these Z electrons fj are distributed 
over j different discrete oscillator frequencies ujj. Thus the oscillator strengths fj must 
fulfil the sum rule ^ • fj = Z. This represents a superposition of more than one of the 
previously discussed resonances. The obtained dispersion relation is: 



e{uj) = H > 9 



£0 m e ^— 4 uj 1 - — uj^ — iwy 



Dielectrics 



Especially for the dispersion of optically transparent substances, as explained, the reso- 
nances do not lie in the spectrum of the electromagnetic waves. Cauchy 15[ replaced the 
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2. Theoretical background 



derived formulae were for this case by an expansion of the resonance behaviour into power 
series with respect to A = 2ticq/uj. The resulting formula is called Cauchy's formula 14 !: 

B C 

y^(X) = l + A+ — + — + ... , 

where A, B, C, . . . are determined from the oscillator frequencies and /,-. Usually in this 
formula only A and B are taken into account. The elegance of this simple model lies in its 
mathematical simplicity, nevertheless, its weak point is that it only describes the dispersion 
relation for regions of normal dispersion. 



To circumvent this deficiency, Sellmeier [87( proposed a description that builds upon 



this approach and includes anomalous dispersion. He expressed the dispersion of dielectrics 
phenomenologically in terms of experimentally determined coefficients Bj for the resonance 
frequencies Uj\ 

Usually the approximation is truncated with m = 4, which is sufficient for a good 
modelling of the dispersion properties of e.g. silica in a wavelength range from the ultra- 
violet, throughout the visible wavelength to the infra-red range (A < 2.3 (j.m). A plot of the 
dispersion for real silica glass (amorphous S1O2), calculated from the Sellmeier equation 
12.91 is presented in fig. 13.11 



Metals 

For a large number of metals, most notably alkaline metals but also for example aluminium 
and silver, the kinetic energy of the electrons is significantly lar ger than the potential energy 
that is caused by the lattice atoms. Thus the Drude model [20| of the free electron gas 
is expected and experienced experimentally to be a very appropriate mode0: Stationary 
atoms, surrounded by a quasi-free plasma of electrons. 

Taking this model into account, the total complex dielectric permittivity of metals can 
be expressed in terms of two separable constituents, the contribution of bound electrons 
Abound and the contribution of the free electron g follows: 



= £bound(w) + £frecM (2.10) 

The contribution from the free electron gas can be calculated analogous to the derivation 
of the dispersion of the bound electron oscillators. The equation of motion of free electrons 

14 A detailed derivation of Cauchy's and Sellmeier's formula, aside from the mentioned original publications, 

can be found in Born and Wolf [l3| . 
15 For a well explained and in depth study of the optical properties of metals it is referred to Dressel and 

Griiner [H], ch. 5, 12. 
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is comparable to that for bound electrons (equation 12. 1.3|) . solely without a harmonic 
potential: 

—eE(r, t) = m c (r + jr) , 

where the damping factor has usually the order of 7 ~ 1 x 10 14 s , which can be intuitively 
interpreted as a relaxation time of the free electron gas (r = I/7 ~ 1 x 10 _14 s). For 
the homogeneous solution of this differential equation the following expression is simply 
obtained: 

1 -it 

7 

which is the description of an exponentially decelerated motion with the decay time r. An 
inhomogeneous solution for a harmonic wave E = Eoe~ lujt represents a periodic motion of 
the electron gas. For N electrons this gives rise for a current density j = Ner. Applying 
the constituent relation j = oE (see p. [8} to j(u>), also the conductivity cr(u>) can be 
obtained: 



r(u>) 



m e (uj 2 + 17a;) 

Ne 2 „ 
m e (7 - luj) 

Ne 2 
m e (7 - iuj) 



E 



For low frequencies w < 7, a can be approximated by its static value of oq = Ne 2 / (m e 7) 
and becomes real. In this case the following equation for the relation of the imaginary part 
of the dielectric permittivity and the conductivity is valid: 



lm(e(uj)) = (JEq/oj 



For high frequencies w>7, which is the case for optical frequencies, a becomes largely 
imaginary (compared to its real part) and the mentioned relation becomes invalid as the 
real and imaginary contribution to e are not separable any more. 

For low enough frequencies the major contribution to the total dielectric permittivity e 
in equation (|2.10[) is caused by the free electron gas in metals. The contribution of the 
bound electrons can then be neglected. From the dispersion relation of a via equation 
(|2.10p and the relation Im(e(k/)) = ueq/uj the total dielectric permittivity of metals in the 
Drude model can be derived: 
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£mctalM = 1 " 2 P , (2.14) 

Ne 2 , . 

w P = , 2.15 



where, analogue to equation (|2.8p for the high frequency limit of the Lorentz oscillator 
model, Wp is the plasma frequency of the free electron gas. For the frequency range below 
that plasma frequency to <C u p , metals retain their metallic, opaque and reflective character. 
For high enough frequencies however, the Drude free electron model's dielectric permittivity 
£f ree is described sufficiently by (|2.8|) . applying the metal plasma frequency. In the very 
high frequency limit eventually lim^^oo £f rcc (u;) = 1. 

For many metals, also those under consideration in the regarded experiment, not only the 
response of the free electrons, described by the Drude model, has to be taken into account 
but also electron excitations between different bands, the so called 'interband transitions] 16 !: 



EmetalM = + Winter M 



Due to their specific electron configuration (the d-band, close to the Fermi surface is 
filled), this effect in the visible range arises especially for noble metals for frequencies 
<jj > uj p . For frequency ranges where these interband-transitions occur the Drude-model is 
not valid anymore. 

Values for the plasma frequencies of typical metals, derived from experimentally deter- 
mined electron densities N [_| are shown in table 12.11 A theoretical treatment of the effect 
is besides the scope of this worlj^l. Usually in these cases experimental data is applied (see 
ch. [3|. 

An analysis of the band structure of gold reveals two interband absorption edges in the 
visible regime, which explain the yellowish colour of the metal. In contrast, the band 
structure of silver does not allow interband absorption in the visible, causing the typical 
neutral colour, which together with the high reflectivity of silver is the reason for its 
widespread application as material for mirrors. 



16 analogue to the 'interband transitions' for metals, for semiconductors 'intraband transitions' occur which 

are not discussed within the scope of this work. 
17 A comparison to theoretical and electron energy loss (EEL) studies can be found infliil]. 
18 Concise derivations can be found in Ashcroft and Mermin [J]. 
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Metal 


Valency 


N in 10 28 m~ 3 


w p in 1 x 10 16 Hz 


A p in nm 


Ag 


1 


5.86 


1.36 


138 


Au 


1 


5.90 


1.37 


138 


Al 


3 


18.1 


2.40 


79 



Table 2.1.: Plasma frequencies of silver, gold and aluminum from [24|. The pl asma frequencies co p 
are calculated with experimental data for the electron density N |l09l | from equation 

Em 

Kramers-Kronig relation 

The real and imaginary parts of the dielectric permittivity e(uj) are in fact not independent. 
The relation between both is described by the Kramers-Kronig relational which shall be 
briefly explained. 

From the frequency dependency of the dielectric permittivity s(uj) the frequency, there- 
fore time dependency enters the correlation of E to D. By a Fourier decomposition of 
D(r,u) into its monochromatic components and a Fourier transformation into the time 



domain, the following relation can be derived [3ll . pp. 381 



D(r,t)=e 



OO X 

E(r,t) + J drG(r)E(r,t-T)\ (2.16) 



V 

oo 



G(r) = h I dw ( £ M- 1 ) e_w ( 2 - 17 ) 



This equation describes a retarded dependency between D and E in space, G being the 
Fourier transformed of rj = e{to) — 1. The causality principle imposes that the response of 
a medium at a time t at a spacial coordinate r can only result from a field that existed 
before t at r. With this premise the lower boundaries of the integrals in equations (|2.160 
are chosen to be and not — oo. 

A corollary of equation (|2.16p is that the dielectric permittivity is an analytical function 
in the complex upper half plane (Ihi(lj) > 0), e(— uj) = e*(uj*). By partial integration of 
equation (|2.17p an asymptotic series is achieved: 



iG(0) G'(0) 

£(U>) — 1 = k h • • • 



19 The Kramers-Kronig-relation was first derived by H. A. Kramers (1927) [4(| and R. L. Kronig (1926) [43 
independently. An in depth explanation can be found in Dressel and Griiner pij 
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Physically in this series the first summand must be 0. Thus, for the high frequency to 
limit Re(e(uj)/eo — 1) decreases with the order of 0(uj~ 2 ) and Im(e(w)/eo — 1) decreases 
with the order of 0(w~ 3 ). 

As the previous corollary is valid, the application Cauchy's integral theorem is allowed 
for relating the real and imaginary parts of e(co) along the real axia 20 !: 



Experimentally the Kramers-Kronig relation predicates that for the determination of the 
dispersion of the complex permittivity e{to) it is only necessary to measure the absorption of 
the medium, which is connected to the imaginary part (Im(e(u;))). The real part (Re(e(o;))), 
can then be calculated. However it is important to point out, that for that objective the 
measured data has to cover a broad spectral range and the sample must not be transparent. 



2.1.4. The vector and scalar wave equations 
The homogeneous Maxwell's equations 

In most optical applications the Maxwell equations can be significantly simplified by the 
assumption of isotropy of a dielectric material or a metal. In both cases the free charge 
density p = oE3 becomes zero. The electric current density in a metal j = oE depends on 
the electric conductivity cr of the material and the electric field E, in the case of a dielectric 
medium j = 0. 

For the following applications additionally a description of homogeneous media is suf- 
ficient. In this case, the dielectric function e as well as the magnetic permeability fj, are 
translationally invariant in space. Under these assumptions Maxwell's equations (12.ip can 

20 Where P represents the Cauchy principal part that is defined for a function f(x) with a singularity at 



DC 





x by: 




21 



Also for a metal with no external field as well as for a metal with a time oscillating field p — 0. 
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be written with divergence free D and B homogeneously: 



V • D 

d 

VxH- e e—E 
at 

d 

V B 







(2.18a) 







(2.18b) 











(2.18c) 
(2.18d) 



The vector wave equations 

Prom equations I2.18bl and I2.18cl the vector wave equations, second-order lineaiF^I partial 
differential equations, for E respectively H are derived, which are also called Helmholtz 



where k is the complex wave number, being the absolute value of the complex wave vector 
k = \k\ of a propagating electromagnetic wave. 

Its imaginary part represents the absorption of the wave in the medium Im(fc(w)) = a/2 
where the intensity of the plane wave decreases as 1(1) = Ioe~ al (Beer-Lambert's absorption 
law). 

In transparent media (Im (e(w)) ~ 0) the wave number can be assumed to be real. From 
equations (|2.19p it is possible to express the dependency of k on the angular frequency uj 
by the dispersion relation: 




V x (V x E) = iujfiV x H = uj 2 e^E V x (V x H) = -iueV x E = uj 2 e^H 
AE + k 2 E = AH + k 2 H = (' 





(2.20) 



The phase velocity of light and its dispersion can then be expressed as: 




IMt c 



(2.21) 




Therefore solutions can be linearly combined to obtain new solutions. 
With the standard notation of the Laplace differential operator AA 
relation [27201 



= V ■ (V ■ A) and the dispersion 
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Plane waves 

The simplest, nontrivial solution of equation (|2.19[) describes propagating 'plane waves'. 

From solving (|2.19[) as well as Maxwell's equations f|2.18[) solutions for the electric and 
magnetic field are obtained. They contain a harmonic (sinusoidal) time dependency (t) 
and a dependency on the spatial coordinate r: 

E(r,t) = £7 e i(fcr ~ a " ) 
H{r,t) = if e i(fcr ~^ ) 

From the divergence equations (|2.18a[ I2.18d[) it can be seen that E • k = and H ■ k = 0. 
Therefore, plane waves, have an orthogonal trihedron of the field vectors and the wave 
vector k _L E _L H _L k, which is the reason to call them also transverse waves Q. 

Particularly for a complex wave vector k = k re + ifcj m the imaginary part fc; m determines 
the amplitude of the field vectors EQe~ kiui ' r respectively Hoe~ kimT . The real part k Te 
determines the phase of the wave tp = k rc ■ r — cot at a certain spatial and temporal 
coordinate. 

Therefore fc re is perpendicular to surfaces of constant phase and fc; m is perpendicular to 
surfaces of constant amplitude of the wave. In vacuum, k re \\ k\ m . The plane wave is then 
called homogeneous, otherwise inhomogeneous. 

The surfaces of constant phase (fe; m = 0) propagate through the space with phase velocity 

v (H2U). 

The scalar wave equations 

The problem of solving the vector wave equations generally, particularly for different coor- 
dinate systems, can be reduced to solving a scalar wave function, which is an easier task. 
Therefore (following ) a vector function M El is constructed from a scalar function ip 
(subsequently called the generating function), and a constant vector c, called pilot vector, 
that can be written aa^l: 

M = V x (c0) V • M = 

AM + k 2 M = V x (c(A^ + fcV)) (2.22) 

24 Due to their importance and as it will be used for the numerical calculations in section [5^4l in appendix lA.2l 
the Poynting vector S is explained that describes the energy transport and direction of electromagnetic 
waves, particularly for plane waves. 

25 Those vector functions subsequently are called vector harmonics, harmonic solutions being the term for 
solutions of the Laplace equation. 

26 applying the identity V ■ (V x A) = 
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M satisfies the vector wave equation if ip solves the scalar wave equation: 

At/; + k 2 ip = (2.23) 

The choice of tp depends on the coordinate system in which the problem is described as 
it has to satisfy the wave equation in that particular case. The vector c then points into 
a particularly exposed direction in the geometric symmetry of the coordinate system. For 
the simplest nontrivial case, for plane waves, the solution is represented by: 



where the wave number, the absolute value of the wave vector k = \k\ is defined for an 
angular frequency to by its dispersion relation in equation (|2.20p . 

M is perpendicular to c. A second vector harmonic N can be constructed, which also 
satisfies the vector wave equation, both vector harmonics fulfil the Maxwell requirements 
(eq.s I2.18bl I2.18c|) for an electromagnetic field: 



V x N = kM V x M = kN (2.24) 

2.2. Planar surface plasmon polaritons 

This work aims for an examination of exciting surface plasmon polaritons on cylindrical 
samples. For a theoretical study of the SPP dispersion relation on the surface of wires, first 
the dispersion relation in the planar case is derived which is then applied to cylindrical 
geometry. 

2.2.1. Background 

The nature of long-range Coulomb interaction in the free electron gas of metals allows the 
excitation of collective vibrations of the electron plasma with respect to the crystal lattice 
by incident electromagnetic waves. These oscillations are called plasmons. The frequency 
of such plasmons in the bulk volume is the plasma frequency of the metal according to 



the Drude model [20J lo p , as previously derived in (12 . 151) . The direction of vibration of the 
electron gas for volume plasmons is longitudinal; the plasma frequency for metals usually 
lies in the range of lo p ~ 1 x 10 16 Hz (see table [2TT]) . 

In presence of a material interface, thus a boundary of the metal, the symmetry is broken 
and new modes arise. This effect was in princip le al ready described mathematically in the 



early studies of Sommerfeld [92| and Zenneck [1121) I . First experimental observation of 



27 These works were already referred to by Debye [17| in his theoretical work on scattering on cylinders that 
are mentioned in section POl 
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dips in the spectra of metal gratings, caused by those 'Sommerfeld's surface wav es' h ad 

13 



m 



already been reported as 'anomalous intensity dips' without explanation by Wood 
1902. These theoretical and experimental works were not combined until Fano [23[ in 1941. 

Shortly after Fano's conclusion, Ritchie [7f| observed unexpectedly low loss in perpendic- 
ular incidence electron diffraction experiments on thin metal foils, followed by several ex- 
tensive analogous experimental studies 7l|, [72J . The results could be explained completely 
by a theoretical derivation of Stern and Ferrell [95|, who called the observed quantized 
excitations 'surface plasmons' (SPf^l. 

After a long time of extended ex perimen ts in the visible wavelength regime especially 
with different prism coupling setups 64 , 100| and extended experiments on the excitation of 
surface plasmons with electron guns, current research on planar surface plasmon polaritons 
focusses on the excitation via surface structuring. This approach is supported by the tech- 
nical progress in microstructuring techniques. Particularly the field of surface plasmonics 
has attracted special interest after Ebbesen et al. 21[ reported on extraordinarily low loss 
in sub-wavelength hole arrays in a metal substrate. Thus also the t heor etical description 



of SPP on structured surfaces recently attracted particular interest 110( | 



There is more than one possible type of surface polaritons. An electromagnetic wave 
that passes a medium induces a polarization (as described in sections 12.1.21 12 . 1 .3[) which 
couples back and therefore by a superposition of incident and emitted field modifies the 
total field. This coupled excitation is called polariton. Surface polaritons (SP) are those 
with the characteristic that the wave is bound to a material surface. 

The SP that this work's aim is to examine are excitations of the electron gas, thus called 
surface plasmon-polaritons (SPP). In general also coupling to phonons could occur, in that 
case leading to surface phonon-polaritons. 



2.2.2. Dispersion of planar surface plasmon polaritons 

In the following, the the dispersion relation of surface plasmon polaritons (SPP) on a 
single, planar interface is derived. An interface is observed between two nonmagnetic 
media 1 and 2 with the wavelength dependent complex dielectric permittivities £i(w) and 
£2(w). The origin of the orthogonal coordinate system (as depicted in fig |2.3p is chosen 
that the coordinate z is perpendicular to the material interface at z = {z < 0: medium 
1, z > 0: medium 2). 

A harmonic plane wave in this system can be described in terms of an electric field vector 
E and a magnetic field vector H. As both materials are taken as free of external charges 

28 Several later re view s and b ooks resume the latest theoretical and experimental developments and the basic 
principles [o9l . Illll . |73| . l8p| . [l06l . loll . IFll. [53. Irj [irj] of surface plasmons (SP) and surface plasmon polaritons 
(SPP) in planar geometry, spheres and cylinders as well as different optical coupling, excitation and 
measurement methods. 
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Figure 2.2.: Schematic of a surface plasmon polariton with propagation vector (3 on a single planar 
metal to dielectric interface, z is the normal vector. The TM field with amplitude \E X \ 
penetrates both adjacent media, decaying exponentially. Schematically the alternating 
SPP charge distribution with correlated fields is indicated. 



p = and currents j = 0, E and H are defined as solutions of the homogeneous Maxwell 
equations (j2.18j) . 

Usually the solution vectors E and H are classified with respect to the surface direction 
into p-polarized light (E \\ x — z — plane) and s-polarized light (E _L x — z — plane For 
waves that travel along the interface, a component of E is necessary that is perpendicular 
to the x — y — plane. Therefore no s-polarized (TE) surface plasmon polaritons exist. Only 
those field vectors are regarded with the additional constraint to be bound to the surface, 
thus the fields must be evanescent for z > and z < 0. The remaining solutions for the 
arbitrarily chosen propagation direction x, can be written for both media i = 1, 2 as: 



E, 



Hi 



( E i>x \ 


( °' 1 



e — K»|«l e — i(/3x— wt) 



g— «i|z|g— i(Px-u>t) 



(2.25) 



(2.26) 



where (3 is introduced as the propagation constant, which is the wave number k of the 
surface plasmon, propagating along the surface. As (3 is defined on the surface, describing 
the SPP on the surface, this wave number must be the same value in both media. From 
the Maxwell equations (12. 181) with (|2.25p . (|2.26p the following relations are obtained: 



p-polarized light generally is light, whose E || plane of incidence (p=parallel), which is also meant by 
the term TM-polarized (transversal magnetic). For s-polarized light E _L plane of incidence, which is 
also known as TE-polarized light. The latter etymologically derives from the German word 'senkrecht' 
for perpendicular. The plane of incidence is always oriented perpendicular to the surface in which the 
k vector of the light lies. 
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Figure 2.3.: Dispersion curves for planar surface plasmon polaritons on a single metal dielectric 
interface, the metal dispersion being described by the Drude model (|2.32[) . The black, 
solid line below the light line indicates the real part of the SPP dispersion Kc((3(u>)) for 
a metal to silica glass interface, whose high frequency asymptote ujsp (|2.33p is plotted 
as a dotted line and whose low frequency asymptote is the silica light line for free wave 
propagation in the dielectric. Analogue, in grey the dispersion and limits for a SPP on 
a metal to air (£2 = 1) interface is shown. For both cases the branches above the light 
lines represent radiative modes. 



I|2.18bll TT i2.1Sb) ,_, . 

iK 1 H lt y = +uje eiE 1}X iK2H 2 , y = -lue^E^x (2.27) 

and the following relation for the propagation constant can be written from (|2.25p . (|2.26l) . 
interpreting k as a z component of the propagation vector: 



UJ 2 



f3 2 - Si^r (2.28) 



c 



Until this point the fields in both media 1 and 2 have been treated separately. However 
on the interface at z = the electromagnetic waves occur a discontinuity in the material 
properties. Due to the law of the conservation of energj@ it can be derived that on material 



30 The proof results in a closed integral perpendicular to the material interface plane as shown e.g. in [96 
ch. 1.3]. 
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boundaries the tangential (x and y) components of E and H must be continuous. Applying 
this condition to equations (|2.27p . the result is the following linear system of equations: 



+ Ho y % = 0, 





which is solved, if the determinant is zero which is the resulting retarded SPP condition: 

— + — = (2.29) 

K\ K 2 

Firstly for a physical interpretation, K\, k 2 are assumed to be real (G M). This equation is 
only fulfilled , as £2(^)1 being the dielectric permittivity of a dielectric is positive and real, 
if £i(u>) of the other material is real and negative. This is true for metals in a frequency 
range where no volume electromagnetic waves can propagate. 

From equation (12.290 the dispersion relation for the propagation constant of planar 



surface plasmon polaritons 76[ is obtained: 



*w(») = «») = r ^? ( T£ ( M < 2 - 30 > 

c y £i{u) + e 2 (uj) 

This dispersion relation can also be written with respect to the free space wave number 
ko = hquj/cq (no = n vac = 1), in terms of an effective refractive index n e fj of the medium 
for SPP: 



y ei(w) +£ 2 (uj) 

For the previously derived dispersion E me tsLi(u) for metals in the Drude model (|2.14p . the 
SPP dispersion relation (|2.30p can be rewritten in terms of the plasma frequency lo p (j2.15|) 
of the material: 





1 ,2 


2u 2 





As depicted in fig l2.3| for a metal (e\{uj)) to dielectric {e 2 {uj)) interface, oj{(3) from equa- 
tion (|2.30p . in the limit of ~ 0) the SPP dispersion approaches the dielectric light line 

— 1/2 

cok£ 2 , which describes the dispersion of an electromagnetic wave, freely propagating in 
medium 2. For higher wave numbers /?, the SPP dispersion function increases monotoni- 
cally but never crosses the light line. 
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Therefore planar SPP can neither radiate light into the dielectric medium nor can they 
be excited by light, incident through the dielectric. 

This conclusion is at least true for perfectly plane interfaces. The excitation of SPP can 
nevertheless be achieved by either surface roughness or gratings or attenuated total reflec- 
tion (ATR) 64], |42[ from inside e.g. a prism whose internally reflecting surface is situated 
near enough to the surface of the metal that light with tunable k component parallel to 
the interface, can be coupled into SPP. 



Although throughout recent literature the terms surface plasmon and surface plasmon 
polariton are used almost interchangeable, in its precise definition a difference shall be 
pointed ouir^l: 

1. For low wave numbers (3 strong coupling between the electromagnetic wave and the 
polarization of the medium occurs, retardation plays an important role, which means 
mathematically that the finiteness of the speed of light cq is taken into account. The 
derived electrodynamical solution (|2.30p then describes a surface plasmon-polariton 
(SPP), that travels along the material interface. 

2. If the propagation constant /3 is much larger than the corresponding plasma frequency 
of the metal (see eq. (|2.15p ) k S> uj p /c, the electromagnetic wave is sufficiently 
mismatched from the induced polarization that limited coupling occurs. In this non- 
retarded limit, the SPP condition (|2.29p reduces to e± + £2 = 0. Non-propagating 
collective vibrations of the electron plasma near the metal interface occur. The 
phenomenon is then called a surface plasmon (SP). 

Strictly speaking, the SP approximation is valid as long as the wave number lies in the 
range of lo s /c <C k <C kp where k-p represents the absolute value of the Fermi wave vector 
(Hsj ] . In fig. 12 .31 it can be seen that in the nonretarded limit for high k the SPP approaches 
the classical, nondispersive SP frequency that with the Drude model from equation (|2.32D 
can be written, as Ritchie [7f|] showed in his pioneering work as: 



wsp = uj p /^l + e 2 (2.33) 

The upper branches in fig. 12.31 represent radiative modes for both material combinations, 
describing electromagnetic waves with frequencies larger than the plasma frequency of the 
metal u > uj p , propagating through the metal. It becomes transparent then. As it can be 
also seen from fig. 12.31 there exists a forbidden range of wave vectors (3 which reaches from 
<^SSP to the plasma frequency uj p . 

In the ansatz for the derivation of the dispersion relation for SPP (|2.25I2.26|) the value k 
was introduced to insert the assumption of exponentially decaying field strengths with \z\ in 

31 This is emphasized and presented in detail by Welford [l06l ]. 
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both media. By applying the dispersion relation (|2.30p to equation (|2.28p . this parameter 
can be expressed as: 



c \ ei+e 2 

This relation allows to determine the so called skin-depth 5{ = 1/k{ which describes the 
e _1 depth of penetration of the electromagnetic fields into the adjacent media (see fig, \2.3\i . 
It is simple to show that the field penetrates deeper into the dielectric than into the metal 
by applying the dielectric permittivities, thus 5\ < 62- For this reason for one metal and 
frequency, the field penetrates deeper into air than into silica glass. 

For a comprehensive description of the propagation of SPP on real metal to dielectric 
surfaces, another effect cannot be omitted. As the dielectric permittivity of real metals 
has a complex part Im(e) 7^ 0, the propagating SPP undergoes attenuation. Regarding 
equation (|2.29p this enters the description in terms of an imaginary part of the wave vector. 
The e _1 decay length of the SPP is then: 

q 1 = 21m (/3(c)) " / £ i^H 



2.3. Surface plasmon polaritons on a wire 
2.3.1. Background and validity 



Schmidt et al. [84[ have recently presented and discussed an approach to convert the disper- 
sion of SPP from the planar case to cylindrical geometry. This model was first introduced 
by Miziumski Hsl, H3j (based upon Pfeiffer et al. 67]) without verification of its validity for 



very high mode orders (m > 40). Also Slater and Frank [90J], pp. 524-537 already proposed 
the approach. 



Liaw and Wu [46J have studied in depth SPP on bent metal dielectric interfaces and 
numerically demonstrated that the dispersion of SPP modes on cylinder like structures 
generally for high enough radii converges to the solution of planar SPP, which in principle 
confirms the validity of the approach. It shall be mentioned that the dispersion relation 
of SPP on lossless circular cylinders is the subject of discussion of several publications 



[62j, [38[. Also the case of metal coated dielectric cylinders is the subject of discussion of 
several recent publications [sH]. 

In addition to the planar SPP solutions, in other geometries e.g. on very small part icles , 



so called localized surface plasmons (LSP) are well known and experimentally proved [111 
p. 142]. They are solutions of the wave equation for the non-retarded SP approximation 
for appropriate boundary conditions and have resonance frequencies that are therefore 
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determined by the size and shape of the particles and the dielectric permittivity of the 
particle material. These non-retarded solutions are valid as long as the size of the particle 
is small compared to the wavelength of the electromagnetic wave (D < A). 

Particularly on the surface of cylinders with very small diameters D<A and very high 
length to diameter aspect ratios I S> A SP excitations arise that share properties of planar 
SPP a nd L SP. D <C A allows non-retarded SP to be excited, I 3> A gives rise to retarded 



modes The result are discrete SP modes on the surface that approach the solution 

of planar SPP (eq. 12 .30|) for high orders. The solutions have both, radiative and non- 
radiative branches of modes for different frequencies, except for the lowest mode which is 
purely nonradiative. 

Thus in the examined case of relatively large radii (compared to A), the assumption to 
apply the proposed model seems valid, at least for frequencies below the plasma frequency 
w p of the metal, as Novotny and Hafner [gl] reasons, a case that is separately discussed in 



detail by Martinos and Economou 54 1. 



However, in this work the derivation of this model for the dispersion and the incident 
angles of excitation for the SPP modes on the surface of a metal cylinder is shown. The 
result is compared to an exact solution of Maxwell's equations for the mode dispersion in 
sec. I 



2.3.2. Model for the dispersion relation 

For a model of the propagation of surface plasmons on a cylinder, the cylinder is treated as 
a closed two dimensional surface on which the SPP propagates as described in section 12.2.21 



Following Schmidt et al. [84J] , the wave vector k^p of a surface plasmon, propagating on 
the surface of a cylinder can be geometrically decomposed (see fig. 12 into a tangential 
component ku, and a longitudinal component (3. f3 = |/3| respectively is the propagation 
constant of the surface plasmon along the axis of the cylinder: 



|fc S pp| 2 = |Av| 2 + |/3| 2 

P = \Aspp-*v ( 2 - 34 ) 

The surface of the cylinder is approximated by a plane tangential to the cylinder surface 
in each point r. Thus the dispersion relation of a planar SPP along an interface of material 
1 inside the cylinder (si = E\) and the surrounding material 2 (e = £2) from eq. (I2.30H is 
applied to eq. (j2.34j) : 
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Figure 2.4.: Schematic of plane waves with the wave vector k- lnc and the angle of incidence a hitting 
a metal cylinder with a diameter D. Resonances occur at angles where (3- mc fits the 
parallel propagation vector f3 which is the axial component of the wave vector of the 
surface plasmon polariton fcspp- The perpendicular component k v lies in a tangen- 
tial plane to the surface of the cylinder on which the surface plasmons form discrete 
spiralling modes. 



The model that is derived bases on the assumption that only integral numbers of field 
nodes can form. This means that 2nR = mX where R = D/2 is the radius of the cylinder 
and m is the mode order, which can be shortened to k = m/R. An additional geometrical 
phase change of n per 180° around the perimeter of the wire is added, which results from 
comparison to an exact derivation of the mode dispersion, thus shifting m to (m — 1): 



(m-1) 



R 



Pm(uj) = ko> 



£i{uj)e {u)) 



\ e*(w) 




(2.35) 



Application of the dispersion relation from equation (|2.20p for a wave vector to the 
propagation constant (3 = luti/cq for each mode j3 m to equation (12.35|) constitutes the 
following relation for the effective refractive index n m of the SPP of mode m, as it is 
defined in eq. (12.31[) : 



n m (uj) 



£i{(jj)e (uj) ( (m - 1)A 



^ £i(u) + £ {u) 



Dn 



Ej(uj)£ (uj) 
\ £i(u) + £ {u) 




2(m - 1) 



D 



(2.36) 
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For reasons of convenience the radius R is replaced by the diameter of the cylinder D. To 
examine the surface plasmons experimentally by side scattering it is necessary to derive the 
incident angles for a matching of the wave vectors of incident light fc; nc in the surrounding 
material 2 to the /cspp vectors of the surface plasmon polariton modes. The light occurs 
with an incident angle a with respect to the cylinder axis. Therefore k is decomposed (as 
shown in fig. I2.4|) : 



cos(a m (cj)j - 



fclncM \y/e (u)\ 
( 



a m (uj) = arccos 



1 



V 



£i{uj)e {ijj) ( c\ I 2(m 



\y/ £ o(u)\ \ ^) + ^) \uj 



D 




(2.37) 



It is important to mention the translation of the polarization condition for the excitation 
of SPP on a cylinder. 

Making the step from the planar case to the model of cylinder-SPPs, "rolling up" the 
planar surface, in principle the basic condition does not change as the model's assumptions 
are not relaxed but tightened by additional boundary conditions. 

Consider a plane wave that approaches with perpendicular incidence (a = 90°) that hits 
a point of the wire, never exactly in the xz-plane (fig. I2.5|) but at least slightly higher or 
lower (y ^ 0). 

The model bases on the assumption that a local approximation of the cylinder surface as 
planar is valid. Thus, by definition a plane wave polarized perpendicular to the xz-plane 
for such a tangential surface is p-polarized (TM). Respectively, if it is polarized parallel to 
the xz-plane, for any tangential plane it is s-polarized (TE). 

In section 12.2.21 it is explained, that SPP on a planar metal-dielectric interface can only 
be TM-modes (p-polarized). These can only be excited by incident light with a matching 
field component in p-polarization. 

This condition for SPP excitation is therefore also valid for the cylinder geometry; p- 
polarized (_L to the xz-plane) incident waves can excite the described spiralling SPP modes, 
s-polarized incident waves (|| to the xz-plane) cannot. 

As soon as the assumption of a = 90° abandoned, the principle does not change. Only 
in the case of polarization || to the xz-plane, the incident electromagnetic field vector Ei nc 
can be decomposed into a radial field component and one that remains s-polarized. 
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2.4. Scattering of light by a small cylinder 




case s 

case p 



Figure 2.5.: Schematic drawing of the coordinate system for the derivation of the scattering theory. 

The incident field is described by (Ei nci H- lnc ), the internal field in the cylinder by 
(Ei, Hi), the scattered field by (E sca , H aca ). The total fields outside the cylinder are a 
superposition of the first and the third. 

The process used to probe the surface plasmons on metallic wires is that of scattering 
of a plane, monochromatic, electromagnetic wave with predefined linear polarization. This 
process can be well described in the formalism of scattering theory. 



2.4.1. Background 

An exact solution of the Maxwell equations for scattering of plane waves on spheres of 
the size of the wavelength was first derived by Gustav Mie and Ludvig Lorenz in 1908, 
eponyms of this solution method, which is also called "Lorenz-Mie scattering". Mie par- 
ticularly obtained a theoretical explanation for the colour effects that were observed on 



light, scattering on wavelength sized spherical particles of gold in a colloidal solution 56]. 



The modern formulation of this theory in terms of Bessel functions was first published by 



Debye [17[ along with a consistent interpretation. 



Unlike the scatte ring solution for the Maxwell equations that was derived by Lord 



Rayleigh (NP 1904) [75J being an approximation for much smaller particle diameters than 
the wavelength, the Mie solution is generally also valid for scattering obstacles whose di- 
mensions are larger than A. The two limits of the Mie description are the Rayleigh theory 
for particle diameters of <C A and classical beam optics for 3> A. 

The Mie ansatz for the solution of the spherical problem can be transferred to particles 
of other shape. Of special interest in this work are amongst those cylinders, having one 
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spatial degree of symmetry less than spheres, which means a slightly more complicated 
mathematical description. 

In the case of an incident plane wave, travelling perpendicular to the cylinder axis (see 
section I2.4.5|) , the problem reduces again and the solution is quite comparable to the Mie 
theory for spheres. The first complete solution of this case for the assumption of a dielectric 
cylinder was already derived by Rayleigh 74 fl, a comparable treatment was presented even 
earlier by Hondros and Debye 30 ] . 

This theoretical description not only was cross checked by a large number of experimen- 
tal and theoretical [97J, |77| applications but it was also experimentally proved explicitly 



e.g. by Abushagur and George Recently Martinos and Economou 54J] have theo- 
retically demonstrated the optical absorption spectrum, applying for this special case of 
perpendicular light incidence on small wires for s- and p-polarization. 

She et al. [89| have extended the description of normal incidence to small (0 <C A) cylin- 



ders that are coated. They assumed the cylinder consisting of dielectric material and the 
coating consisting of a metal (or a metamaterial) and vice versa, neglecting dissipative 
effects. The theory is validated by a comparison to an analytic description accounting for 
retardation of the surface plasmon modes, derived by Schroter and Dereux [85j for dielec- 
tric cylinders (0 = 1 |j.m) that are surface coated by a thin layer of metal. 



The first analytical solution of scattering on a dielectric cylinder for oblique incid ence 
following the Mie ansatz was in some abstract formulation presented by Wait 1031 ] who 
also compared this plain wave solution to incident spherical waves [1021 ] . Kerker et al. 



371 ] , Kerker [36( subsequently treated the problem more extensively. In the case of oblique 
incidence the scattering problem cannot be simplified significantly, thus the theoretical 
description becomes more extensive. 



Kuik et al. 44] lately have shown by applying the T-matrix method 33 ! that the scattering 
solution for cylinders of finite longitudinal extension for perpendicular and oblique incidence 
does not converge to the exact solution for infinite cylinders as the aspect ratio (ratio of 
length over diameter) inc reas es, an examination that was with exact solutions also done by 



Wang and van de Hulst 104]. Lock [48] has extended the Mie type solutions for incident 



plane waves on infinitely long cylinders to non plane wave real beams. 

Nevertheless the derivation for infinitely long cylinders and plane waves that is presented 
is assumed to show good agreement with the conditions of the experiment, as the extension 
of the experimental beam, that can be well treated as a plane wave, is large compared to 
the lateral extension of the wires but small compared to the longitudinal extension of the 



Which is well presented in [9611 . A complete historical review about the founders of the modern scattering 
theory can be found in [351 ] . 
33 The T-matrix method (T: transmission matrix) is an alternative numerical method, also called extended 
boundary condition method', that allows for exact scattering computations f or pa rticles with arbitrary 
shapes and optical properties. For details about this method it is referred to [105| ]. 



2.4. Scattering of light by a small cylinder 



31 



wires. Therefore particularly effects of scattering on the wire's edges that were considered 
in the mentioned publication can be omitted in the examined case. 

The following analytical derivation of the scattering amplitudes and field distributions 
from scattering by an infinite right circular cylinder follow the derivations of the theory 



of scattering on particles of di fferent shapes by Bohren and Huffman [12|] and the earlier 



treatise of Van De Hulst 10l| . Both concentrate mainly on the more common case of 
scattering on spheres, thus their in depth discussion of the spherical geometry is in the 
following transferred to the cylindrical geometry. 

The results of a subsequent computation (presented in sec. I5.3|) of the theoretical scat- 
tering amplitudes should be substantially comparable to the experimental results of this 
work as well to the scattering by metallic nanowires embedded in a dielectric matrix as 
to the scattering on coated cylinders. The latter with the assumption of a high ratio of 
the coating thickness to the skin depth so that it can be treated as bulk material for the 
plasmonic effects. 



2.4.2. Cylindrical coordinate system 

The scattering problem is solved in the intrinsic coordinate system of a cylinder, r = 
(r,(p,z), the system's intrinsic symmetry is exploited, thus the complexity of the problem 
reduces and an analytical solution becomes possible. As drawn in fig. 12.51 the cylinder axis 
defines the coordinate axis z, the extension of the cylinder is defined by its radius R and 
its length is taken as infinite. 

The incident light is assumed to be a plane wave whose fcjnr, I lectin vector lies, in cartesian 
coordinates, in the x-y plane, enclosing an incident angle of a r with z. The observed 
portion of the scattered field's propagation direction (fc S ca||e S ca) lies in a plane with the z 
axis that encloses the azimuthal angle ip with the incident x-y plane. 



2.4.3. Field expansion in vector cylindrical harmonics 

The starting point of this derivation is the scalar wave equation (I2.23P that has to be 
written in cylindrical coordinates: 



13/ dip\ 1 d 2 ip d 2 ip 
r dr \ dr J r 2 dip 2 dz 2 



The definition of the incident angle a follows the convention of the experimental chapter^ Note that this 
is in accordance with the convention of Boh ren and Huffman flit ] and most later publications referring 
to this. Different to that Van De Hulst [lOll ] along with nearly all publications earlier than Bohren and 
Huffman [3 defines -n/2 <a< +n/2. 
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In these coordinates the generating function ip is calculated using a separation ansatz to 
be: 



Mr,^z) = Z n (p)e in ^e ihz p = r V / k 2 -h 2 (n = 0,±l,...) (2.38) 

Where h and n are separation constants, h is chosen to fulfil the boundary conditions. 
In this case for a plane wave with an incident angle a, h = —k Q cosa. Z n (p) are Bessel 
functions of integral order n. Their arguments are scaled, dimensionless radii p. 

In the case of cylindrical coordinates ip generates the cylindrical vector harmonics M 
and N as defined in equations f|2.22[> respectively (|2.24jl to bell: 



( m\Z n {p) \ 

-Z'n(p) 





N n = ^k 2 - h 2 \e i( w +h ^ | -hn\Z n {p) 

JW=Wz n { 9 ) 



\hZ' n {p) 

1 - 



r,(p,z 



These vector cylindrical harmonics are orthogonal, forming a basis. M n ,N n are not nor- 
malized pair wise but over the summation of all orders n. 

The electric and magnetic fields involved in this scattering process are divided into 
three domains (see the right part of fig. l2.5p . Without any scattering object the incident 
electromagnetic field would simply pass the space undisturbed. The fields describing this 
wave throughout the whole space with e Q are called E- mc and i?i nc . 

As soon as the material of a cylinder, placed in that space is described by an e\ ^ e Q 
the electromagnetic field inside the cylinder changes and is described by only E{ and H\. 
From the surface of the cylinder a third electromagnetic wave is emitted, which represents 
the scattered fields E scai and H sca . 

To calculate the total field components in the space outside the cylinder, a superposition 
of E Q = E inc + .E sca respectively H Q = H inc + H sca is necessary. 

All electric and magnetic fields are to be expanded in the vector cylindrical harmonics as 
in this representation the boundary conditions are easiest to apply. Two different cases for 
perpendicular polarizations are solved separately, linear combinations of which represent 
all possible polarizations of the incident plane wave. The time dependency of the fields 



35 If indicated as in this equation, vectors are referring to the cylindrical coordinates, their entries corre- 
sponding to r, tp, z. 
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can be expressed as e lujt and will hereafter be omitted as the whole problem can be solved 
statically. For s-polarized incident light (case s), thus E inc \\ x — z — plane the expansion 
results in: 



E inc = E e ik ° r = E 



V 



- cos a 


sin a 



e ifc r = g EnN {i) ■ Eq = Eo 



ko sin a 



Expressing the expansions for each of the fields in terms of En the following relations are 
defined: 



Ei 



E K 



E 

=— oo 
oo 

E E n (f n N$+g n MW 

=— oo 
oo 



-iife 



j/ inc = — e 



-iJfe 



(2.39a) 



n=— oo 
oo 



H s 



ik 



n=— oo 



(3) 



(2.39b) 



For the perpendicular, the so called p-polarization of the incident light (case p, E\ _L 
x — z — plane) the expansion for the incident field changes to: 



oo 

i E 

71= — OO 

2.4.4. Boundary conditions and analytical solution 

From a macroscopic point of view on the boundary between the two materials e± = e; (in 
this case inside the cylinder) to £2 = s Q (in this case outside the cylinder) the electromag- 
netic waves occur a discontinuity in the material properties. For the further discussion an 
abbreviation is defined as E\/e =: m 2 . 

This occurs on the surface of the cylinder at r = R. Due to the law of the conversation of 
energy^! t can be derived that on material boundaries the tangential (tp and z) components 
of E and H must be continuous. 

36 analogue to the derivation of the dispersion of planar SPP on page 1231 Equations 12.401 base upon the 
assumptions of finite current density and finite conductivity which are both true in physically real cases. 
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e r x (E -Ei) = 



e r x (H Q - Hi) = 



(2.40) 



This leads to a linear system of 4 equations that can for case p be written in matrix 
notation as: 



fhnj;J n (r)) r]J' n {ri) 
7]hJ' n (rf) hnJ n (rj) 

rfJnirj) 











ffn\ 




(hn±J n (0\ 


} ^ 






9n 






} 







b n 




rJ 2 U0 


} E, 



View 



(■^inc j-f^inc) 



(2.41) 



For case s, as E- mc is defined different (see equation I2.4.3P the right side of equation 12.411 
changes to: 



-ihnJ n {£) 




11,!, 

E z 



In these equations J n are Bessel functions of first kind and integer order n, J' n 
\{Jn-\{%) — Jn+i{x)) are their first derivative;!^! [94]. Respectively H^P 



J n + \Y n are 

the Bessel functions of third kind (also called Hankel functions of first kind) whose first 
derivatives are H^ l \ n = ^{H^_^{x) — H^^x)). It can be shown that Hankel functions 
of first type H^ 1 describe outgoing and Hankel functions of second type Hn^ = J n — iYn 
describe incoming waves. Thus, it is according that in this case Hn is the choice for the 
description of the scattered wave. 



For reasons of clarity the following two abbreviations are introduced for the limits of p 
(see equation I2.38[) on the cylinder surface from inside the cylinder p x and from outside the 
cylinder p Q : 



37 Useful information about the Bessel functions, the derivatives and the Bessel differential equation are 
gathered in appendix lA.il 
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v 9Jk?- h 2 R = k Q R\/m? - cos 2 a = Pi (r = R-5) (2.42) 

f IP y/kl -h 2 R = k Q R^l - cos 2 q =A: OJ Rsina = /3o (r = J R + 5) (2.43) 

By analytically solving the equation system (|2.4ip relations for the coefficients a n , b n , 
f n and g n are achieved. From these with equation (|2.39|) the electric and magnetic field 
vectors inside the cylinder (2%, Hi) and the scattered field vectors outside the cylinder 
(E sca , H sca ) can be evaluated. 

The solution vector (a n ,b n , f n ,g n ) includes rather extensive expressions. After some 
simplifications and introducing abbreviation expressions A n ,B n ,C n ,D n ,V n ,W n the following 
standard notation^! can be achieved for the two perpendicular polarizations in case p and 
case s for the scattered light outside the cylinder: 



a 



_ C n V n B n D n 
" W n V n + iDl 
A n V n — \C n D n 



W n V n + ID 2 



W n B n + iD n C n 
W n V n + iDl 

. C n Wn + A n D n 

1 W n V n + iD* 



(2.44a) 
(2.44b) 



In these equations against a number of useful abbreviations is introduced, shortening 
the expressions significantly: 



B n = £ (m 2 ZJn(v)UO ~ ^nivKiO 
C n = ncos(Q)7 ? J n (r ? )J„(e)(e 2 /^ " 1) 



The scattering coefficients show symmetries that are caused by the symmetries of the 
Bessel functions J_ n , = (—l) n J n and Y- n = (— ^) n Y n , which was already demonstrated by 
Kerker et al. [371 ] : 



38 This notation can be found in Bohren and Huffman slightly modified in Van De Hulst [ 10 it ] and 
certain recent papers (e.g. [H3. 
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O'—ns — ^ns ^— ns — ^ns <^0s — 

^— ns — 0"np b— ns — b n p &0p — 

<^ns — ^np 

2.4.5. Perpendicular incidence 

The special case of perpendicular incidence (a = 7t/2) is of particular interest as equa- 
tions (|2.44al I2.44bp can be simplified significantly as £ = k Q Rsm(a) = k Q R and n = 
k R\fm 2 — cos 2 (a) = k Q Rm. Thus, the scattering coefficients can be expressed explicitely 
as: 



«np± = 

J n {mkR)J' n {kR) - kRJ' n (mkR)J n {kR) 

bnpL — 7TT7 777 

J n (mkR)H { n> (kR) - mJ> n (mkR)H { n> \kR) 
mJ' n {kR)J n {mkR) - J n {kR)J' n {mkR) 

a ns±. — 7TT7 77T 

mJ n {mkR)H { n> (kR) - J^(mkR)H [ n' (kR) 
b ns _L = 



As the complexity of these simplified equations is much lower than that for arbitrary 
incident angle a the case of perpendicular incidence has attracted special attention in 
literature and was discussed e.g. by Martinos and Economou 54] H 



2.4.6. Scattering efficiencies 

The scattering process cannot only be described in terms of the scattered field distribution. 
The natural measurable scalar quantity of electromagnetic field distributions is the intensity 
I(r) defining the total electromagnetic energy flux per area (and is sometimes also called 
irradiance). The intensity mathematically is derived from the nonobservable Poynting 
vector S of an electromagnetic wave. 



J(r) = |5(r)| S(r) = 1/2 Re (E(r) x H(r)*) 

The Poynting vector includes additional information of the direction of energy fluj^l. 
In general S is a complex vector, whose real and imaginary parts represent the physical 

39 In this publication the coordinate system thus also the scattering coefficients are defined different from 

the demonstrated derivation. 

40 For details it is referred to appendix I A. 21 
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energy flux respectively an analogon to the idle power. In this thesis the Poynting vector 
is defined as the real part only. From S it is possible to derive a more general description 
of the scattering efficiency. 

In precise terminology the scattering process is divided into two contributing effects, 
together forming the total extinction, which is also called attenuation: 

extinction = scattering + absorption 

Both processes, scattering and absorption withdraw energy from the incident electro- 
magnetic wave thus reducing the intensity of the transmitted wave. Scattering describes 
a deflexion of electromagnetic waves, thus a change in the energy flux distribution. Ab- 
sorption comprises the dissipative processes involved in the scattering, including also the 
absorption that is connected to plasmon excitation and the predominantly ohmic losses on 
the SPP. 

The electromagnetic energy flux through a surface A is defined by the surface integral 
over the vector component of S perpendicular to the surface (S ■ e±). The definition of a 
surrounding, closed sphere with surface A' and radius R' > R surrounding the scattering 
cylinder leads to an integral equation for the net energy flux to (W > 0) respectively 
from the sphere (W < 0), the first being omitted in our application of scattering of an 
electromagnetic wave from a source outside the sphere: 

W = - J dAS ■ e r (2.45) 

A' 

For the described system the law of the conservation of energy must be fulfilled. Thus 
the sum of energy flux through the surface A 1 for all processes participating in the total 
scattering process must balance out to 0, generally including absorption. 

Then the energy that is absorbed by the scattering obstacle can be expressed as Wabs = 
Wi nc — Wgca+Wext- For a nonabsorbing surrounding medium (always valid as Im(e(Si02)) = 
and Im(e(vac.)) = 0), without any scattering obstacle the incident waves passed the 
volume element that is restricted by the surface A' without loss of energy. 
As this assumption is true, eq. (I2.45P can be simplified to the following expression for A 1 : 

W cxt = Wabs + W SCSL 

These energy flux values do not yet only depend on the scattering configuration but also 
on the incident irradiation /. Therefore a new variable is defined by dividing this linear 
dependency out, leading to the so called cross sections [C n ] = 1 m 2 : 

Wext 
I'mc 

Cabs C sca 



C ex t 
C ex t 



38 



2. Theoretical background 



From these scattering cross sections of the unit of area a unitless value: The scattering 
efficiency Q n which is defined by dividing the projection of the geometrical cross section 
G (that can be calculated for each particle shape) on a plane perpendicular to the incident 
beam C n as Q n = C n /G. The total scattering absorption efficiency is then given by: 

Qabs — Qcxt Qsca, • (2.46) 

From a strict point of view the transfer of this theory to the case of an infinitely extended 
cylinder is not as obvious to be valid. It is well defined for spatially confined particles that 
can be simply encased by a closed plane. But as the treated cylinder is infinitely long, its 
scattering C sca and absorption cross section C a b s are of course infinite. 

However, the unitless definition of Q n (eq. 12.46]) allows a differential description of the 
scattering and absorption efficiency even for infinitely extended cylinders and infinitely 
extended plane waves. 

Expressions for the absorption efficiency of spherical and cylindrical particles were also 



calculated from the Rayleigh theory. Luk'yanchuk and Ternovsky 49J presents the Rayleigh 
scattering efficiency of an infinitely long cylinder as: 



n 2 



m 2 — 1 



m 2 + 1 



e 1/2 
A 



However this description does only take into account dissipative effects (Im(ei) ^ 0) and 



not radiative damping that is inherently included in the Mie scattering description [49]. 
Additionally the Rayleigh solution, as mentioned, is only applicable for small scattering 
samples (A 3> D). 

The derived scattering theory can also be used to calculate the scattering efficiencies Q n , 
which do not show the weak points of the Rayleigh approximation. Therefore usually an 
amplitude scattering matrix T is derived that includes all effects and connects the incident 
field components (E- mc ,H- mc ) in cylindrical coordinates to corresponding scattered field 
components (E- mc , H- mc ). 

In the case of a cylinder, the surface of integration of the Poynting vector is defined 
as a concentric cylinder with length L', surface A' and radius R 1 [l3]. Integrated are the 
components of the Poynting vector perpendicular to this surface, the radial components 
(S) r : 



» o/r/ f 2U , ../« x w m o/r/ 



W sca ^ R L / d^(S sca ) r W ext R L / d^(5 cxt ) r 
Jo Jo 

S sca = i/2Re (E sca x H* sca ) S cxt = i/ 2 Re {E inc x H* SCSL + £ sca x H* nc ) 
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Introducing the expansions of E sca _ and E- mc (eq. I2.39al I2.39b|) into these equations and 
integrating, the scattering and extinction efficiencies of the infinite circular cylinder are 
obtained: 



Q ^> s = 2RUJ- = k^R ^ s0 ' 2 + 2 g ( IM2 + 

Qsca, p = 7—7^ I | fl po| 2 + 2 ^ ( | fl pn| 2 + |& 



|2 

sn 



00 
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pn 

n=l 



Qcxt, p = ]^R Re [ a P° + 2 Yl 



00 
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Figure 2.6.: Spatial plots of the first orders of the vector cylindrical harmonics. The first row dis- 
plays the amplitudes of the functions Re(Mo) 2 and Re(iVo) 2 in two of the perpendicular 
main planes, the second row for Re(Mi) 2 and Re(iVi) 2 , respectively the third row for 
Re(M2) 2 and Re(A r 2) 2 . All calculations are performed for a wavelength of A = 632.8 nm 
and free space i.e. e = 1. The spatial coordinates x, y, z are displayed in |j.m, the color 
scales arbitrarily normalized. 
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2.5. Asymptotic scattered field and intensity pattern 

In the experiment the observation point is comparably far away from the scattering centre. 
Introducing far field approximations, also the form of the observed intensity pattern can 
be derived 41 !. 

For large distances from the sample, as fco-Rsin(a) ^> 1, the scattered field can be 
asymptotically written as: 



TP -_ R nP ~ i7r / 4 / ik (rsm(a)-zcos(g)) _ 

■ c 'sca,s — - Cy t; \ tj- n • f \ 

y 7tAo-Rsm(a) 
y^(-l)"e my (a sn e^ + a sri (cos(a)e r + sin(a)e z )) 



TP -_ Rnp-" 1 / 4 , I lfcp (r sin(a) cos (a) ) 

-E'sca.p — -E'Oe 4 / _ . , . e ^ 

y nK Rsm(a) 

^(-l) n e my {-a pn e v - 6 pn (cos(a)e r + sin(a)ej) 



The wave fronts of the scattered wave (planes of constant phase) in cartesian coordinates 
are found to satisfy the relation: 



f(x, y,z) = r sin(a) — z cos(a) 

This relation describes cones of half angle a. Therefore the scattering pattern that 
is expected from the experiment is the projection of cones. Concisely for a = 90°, the 
cone is expected to reduce to a cylinder. In a plane perpendicular to z at a = 0, for 
angles a < 45°, the projection of the calculated intensity forms ellipses with decreasing 
eccentricity, approaching a = 0°. For a = 45° the pattern is a parabola and for further 
increased incident angles, hyperbolas are observed. 



Details about the derivation can be found in Bohren and Huffman [13 ] . 
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Chapter 




Sample fabrication 

The following experiments were performed by the using three different materials: Silica 
(fused silica glass, S1O2), gold (Au) and silver (Ag). For the computation of the dispersion 
of the surface plasmons (sec. l2.3.2|) as well as for the computation of scattering field dis- 
tributions and excitation angles (sec. l2.4p the knowledge of the material dispersion e{oS) of 
these materials is essential. 

Of additional interest are the physical properties regarding the fabrication process of the 
samples. Thus, an overview over the relevant information of all three materials is presented. 
With this information background, the techniques that are applied for the fabrication of 
two different types of samples are summarized: 

1. Free standing wires, consisting of a silica fiber substrate with a diameter that is varied 
between 10 |xm < D < 40 [im, coated with a comparably thick layer of sputtered 
metal. 

2. Thin wires with a diameter 0.4 \xm < D < 1 [im which are fabricated by drawing 
a specifically designed fiber, analogue to the fabrication of photonic crystal fibers 
(PCF) with a single hole or a hole array. The metal that is subsequently filled into 
its holes forms arrays of wires inside the remaining silica glass matrix. 



1 As silica is an amorphous solid, no discrete melting point exists. Provided is therefore the softening point 
at 1730 °C, additional parameters which are important for tapering is the annealing point at T = 1075 °C 
and the strain point at T = 1180 °C, which are all defined by different grades of material strength (see 
for more details [281]'). 
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Material Processing properties Optical properties 

T m (°C) T b (°C) a oxp Type Model e(632.8nm) 

Ag 961 2162 18 x 10" 6 n.m. fee Drude -16.1396 + 1.0814i 

Au 1063 2856 14 x 1(T 6 n.m. fee experimental -9.56 + O.Sli 

Si0 2 1730 1 . 2230 5.5 x 10~ 7 amorphous Sellmeier 2.123 



Table 3.1.: Overview of the processing (melting T m (in °C) and boiling point Tb (in °C), thermal 
expansion coefficient a CX p in 1 x 10~ 6 K _1 at room temperature) and optical properties: 
Material type, dispersion model, dielectric permittivity e(A = 632.8 nm) of the noble 
metals silver (Ag) and gold (Au) and of the applied type of fused silica glass (Si02) with 
a pureness of 99.99 %. 



3.1. Material properties 
Silica 

For the sample fabrication pure silica glass is usec@. The amorphous material is fabricated 
by fusing crystalline silica (SiC>2). Silica glass is an insulator (a ps 0) and transparent in the 
visible wavelength regime (see section 12.1.31 p.ITT]). As the SiC>2 molecules exhibit electronic 
vibration resonances in the ultraviolet, the absorption in the UV increases for A < 200 nm. 
Vibrational absorption makes silica intransparent in the infrared for A > 2000 nm. 

The dispersion of silica glass is well described by the Sellmeier equation f)2 ,9f) . which 
practically is often written in terms of dimensionless coefficients Aj and coefficients Aj of 
the dimension of the wavelength as follows: 



e(A) = 1 + 



^4iA 2 
A 2 



A? 



+ 



^ 2 A 2 
A 2 



A 3 



+ 



A 3 X 2 
A 2 



\2 
A 3 



(3-1) 



The Sellmeier coefficients always have to be experimentally determined. In this case 
data for bulk fused silica from Agrawal 0| is taken, who fitted the Sellmeier equation 13.11 
into a measured dispersion curve from Sotobayashi and Kitayama [93J : 
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0.696166 3 


0.407426 


0.8974794 


ujj in m 1 


2.753 703 x 10 16 


1.620465 x 10 16 


1.903416 x 10 14 


Aj in m 


0.068 043 x 10~ 6 


0.116 2414 x 10~ 6 


9.896161 x 10" 6 



As depicted in fig. 13.11 the dielectric permittivity of silica glass in the visible range is 
2 Hilgenberg 0620 with a purity of 99.99 %. Its physical properties from [2|| are included in table [5711 
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approximately e ~ 2.123 (re ~ 1.457), which decreases by 1% in towards the region of 
A w 1.5 iim. 



Gold 

As mentioned in section T2.1.3I (p. fl~2l) . the Drude model is not applicable for gold in the 
visible spectrum, as this metal exhibits strong interband-transitions in this wavelength 
range. Therefore for the computations commonly used experimental data from Lynch 
and Hunter [50( |. pp. 286-295 is interpolated for a wavelength range of 200 nm < A < 
1000 nm. The experimental data in the wavelength range of interest is obtained by Theye 
[9^], measuring the reflectance and transmittance of thin (10 — 25 nm), semi-transparent 
evaporated Au films on superpolished fused silica substrates. 

Recent results of density functional calculations and REELsH measurements on the di- 
electric proper ties of silver (Ag) and gold (Au) in a range from infrared to ultraviolet by 



Werner et al. 1071 ] indicate that the widely used experimental dispersion dataset from 



Lynch and Hunter 50J] deviate from new values. This deviation increases for the high en- 



ergy range. For the applied lower energy range (A = 633 nm : 1.96 eV) the results from [50] 
do not deviate significantly, at least nor new, better experimental neither numeric results 
have been reported. Thus the used experimental values are still the best available for the 
desired wavelength range. 



Silver 

The optical properties of silver are well described by the Drude model of the free electron 
gas as described in section 12.1.31 (p. [T2]) . 

Different to gold, for silver interband-transitions play a negligible role in the visible 
frequency range. The interband absorption edge of silver is in the UV at around A = 310 nm 
(u = 6.08 x 10 15 Hz). 

The complex dielectric permittivity £Ag(w) can therefore be written using the dispersion 
relation (|2.14p . separated into its real and imaginary part and converted to a convenient 
form: 



Re(£ A gM) = A 1 - , (3.2) 

Im( £Ag H) = B 3 + B 2 B n2 » ( 3 - 3 ) 

where the coefficients for silver, calculated from the plasma frequency of silver lo p (as listed 
in table EU p. CED are Ai = 3.9000, A 2 = 1.3325 x 10 16 , ^3 = 6.992 407 x 10 6 , B x = 



3 REELS: Reflection electron energy- loss spectroscopy 
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Figure 3.1.: Material dispersion of silica (amorphous Si02), calculated from the Scllmcicr equation 
(|3.ip . The imaginary part of the dielectric permittivity of silica is negligible, its real 
part is displayed. 
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Figure 3.2.: Material dispersion of gold, interpolated from experimental data [50(. For Au in this 
wavelength region intcrband transitions arc not negligible and the Drude model (see 
section B.l. 3 li is not valid, as is clearly noticeable by the increased absorption (Im(e(A))) 
for A < 600 nm, compared to the dispersion of silver as calculated by the Drude model 
and depicted in figure I3.3I 

4i 1 1 1 1 1 1 1 1 1 20 




Figure 3.3.: Material dispersion of silver (Ag), computed from equation (|3.2p . (|3.3|) for the Drude 
free electron gas model (sec. I2.1.3)l . 



48 



3. Sample fabrication 



1.0907 x 10 17 , B 2 = 1.46511 x 10 12 , B 3 = 0.420 62. A plot of the dispersion, computed 
from these values is shown in fig. 13.31 



3.2. Free-standing wires 
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3.2. Free-standing wires 

For the investigation of SPP on wires with a diameter 10 iim < D < 40 |J.m a coating 
technique was applied to fabricate samples. This method involves tapering of circular 
silica fiber and sputtering with a layer of metal with a thickness large enough [d > 0.5 \im) 
to treat them as virtually bulk metal wires. 

The advantage of this method lies in its flexibility to access a wide range of diameters 
precisely with a high uniformity of the diameter and shape along the longitudinal direction. 
The technique of heat tapering of silica glass fibers makes it is possible to fabricate wires 
with very high aspect ratios. At least sample lengths of 0.4 cm are necessary as the wire 
diameter on the one hand must be much smaller than the probe beam diameter to fulfil the 
condition of plane waves (p. l30l) . On the other hand for a variation of the incident angle, 
the effective projected beam diameter into the longitudinal direction of the wire increases 
with low incident angles (compare to fig. l4.ip . making a sufficient length of the sample 
necessary. 

The film preparation of these fibers is performed using the magnetron sputter deposition 
technique, allowing well reproducible results with a high surface uniformity and low surface 
roughness. 



3.2.1. Tapering 

As starting material for the core substrate silica glass fibers with diameters of D" = 125 |xm 
or D" = 200 (xm are taken. These fibers are stripped of their protective coating and cleaned 
accurately with a cascade of acetone (CH3COCH3), isopropanol (CaHgO) and ethanol 
(C2H6O) to achieve a clean surface without any residuals. 

In the next step (fig. l3.4l a) the diameter of the fiber is tapered down by the application of a 
butane-oxygen flame that is continuously moved along the fiber with a computer controlled, 
well defined speed and increasing movement distance while a controlled tension is applied 
between both ends of the fiber. The temperature of the fiber for this process is kept below 
the temperature where the silica's viscosity decreases to far, thus its circular cylindrical 



shape is well maintained. This technique is commonly used [34| and well described e.g. by 
Birks and Li dj. 

The fiber diameter D" for a certain length I is uniformly tapered down to a smaller 
diameter D' with a fast transition region to the original fiber diameter on both sides. The 
fluctuations of D' over the tapered length are measured with an optical microscope to 
be AD' < 0.01 \im. For the sample fabrication usually a length of I ~ 2 cm from the 
center region of a uniformly tapered length of 4 cm is made at a final fiber diameter of 
10 (im < D' < 40 (xm. 



4 Drawing c and the magnetron sketches are reproduced following [K 
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Figure 3.4.: Schematic drawing of the tapering and sputtering process for sample fabrication. Part 
(a) displays the flame tapering process, reducing the diameter of a fiber sample from D" 
down to D 1 . This sample is then sputtered with a virtually monoatomic contact layer 
of Cr and a thick layer of Ag or Au by magnetron sputter coating. The final, virtually 
bulk metal sample has a diameter of D. Inset (c) illustrates the dependency of the 
corning size and shape on the sputtering temperature T with respect to the melting 
temperature T m of the metalQ 



3.2.2. Sputtering 

The tapered fibers with diameters D' are vertically mounted on the rotating turntable 
inside the vacuum chamber of a sputtering system_|, which is equipped with more than 
three magnetron sputtering devices. 

Between the silica glass substrate and the deposited silver or gold layer always a very 
thin contact layer of Cr is sputtered to achieve a better adhesion between both materialqj. 

Inside the chamber a vacuum of (1 x 10 -6 bar) is established with a turbomolecular pump 
system. A low pressure argon atmosphere (1.3 x 10~ 6 bar) is then conditioned in which a 
gas discharge ignites a continuous plasma. Argon ions from this plasma are accelerated 
by the magnetic field of the magnetron devices and impact into the sputtering target (Cr, 
Ag, Au) where neutral atoms of the sputtering target are knocked out. These single atoms 
travel with the starting energy towards the sample surface where they assemble to form a 
metallic layer on the glass substrate. 

5 AJA OTC Orion series uhv, also used was a smaller single magnetron sputtering device Emtech K575X. 
6 This work shall only give a short summary of the applied magnetron technology. For a more detailed 
treatement the reader is referred to the very good introductions [H^] and the recent lecture script [Io| . 
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Figure 3.5.: Scanning electron micrograph and FIB cut into the surface of a sample with diameter 
22.64 um and a thickness of the sputtered gold layer of 670 nm. Insets 1 to 3 show SEM 
micrographs of the surface with different magnifications. Insets 4 and 5 depict the 
result of a focused ion beam cut perpendicularly into that surface. During the cutting 
process a thin layer of carbon (C) was deposited on top of the gold (Au). Clearly the 
structure of the gold layer is visible. The nm range thin contact layer of Cr below the 
Au layer is not visible. 



The lower the pressure inside the chamber, the higher is the purity of the sputtered 
material and the higher is the energy of the impinging sputtering molecules on the sample. 
On the other hand, as the gas pressure decreases below a limit, the plasma extinguishes 
and interrupts the sputtering process. Sputtering deposition offers generally the advantage 
of comparable high surface smoothness. 

To achieve a minimum surface roughness (A), thus grain size, a precise adjustment of the 
sample temperature with respect to the melting point of the sputtered material (T/T m ) is 
essential as is illustrated by fig. 13.41 inset (c) . Too high temperatures support the formation 
of large grained material spheres on the substrate surface (D). A compromise was chosen, 
also taking into account that the very thin silica substrate fibers are very sensitive to high 
temperatures in the geometry as they are perpendicularly mounted in the sputtering cham- 
ber. 
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To monitor the surface quality, usually examination with an optical microscopd_| was 
sufficient. For adjusting the sputter deposition paramers, using a focussed ion beam (FIB) 
and SEM systenjj, a series of micrographs of the obtained samples was taken as shown in 
fig. 13.51 In pictures (1-3) micrographs of the surface of a sample with an outer diameter of 
22.64 urn are shown. The sample exhibits a rather small corning size of < 100 nm. 

For the examination of a cross section of the sputtered layer a thin layer of carbon 
is first deposited on the sputtered layer of gold (Au), allowing a perpendicular fib cut 
with accelerated gallium ions into the layer, whilst inevitably the removed matter deposits 
around the cutting place. Insets (4,5) show the resulting view of a cross section of the 
sputtered Au layer. The layer thickness is, as calculated from the assumed sputtering 
rates, 670 nm. A distinct structure inside the layer is visible with thin longitudinal hole 
sizes up to half the layer thickness. Successive cutting and SEM observation of the layer 
disclosed that these structures do not show a large lateral extension. 

Nevertheless it is assumed that the observed structure degrades the SPP excitation 
properties of the sample. A further reduction of the hole structure of the sample layer could 
therefore decrease the impurity level of the metal, hence further increase the conductivity 
and support the SPP excitation and guidance properties. 



3.3. Wires in fiber silica matrix 



To fabricate metal wires of very low diameters and very high aspect ratios a technique is 



applied that was recently published by Schmidt et al. [84|. Silica glass fibers are drawn 
that provide a hole structure with longitudinal extension in a specifically defined order and 
size. These holes are then filled with metals, which form wires of very small diameters 
(D > 400 nm) and extraordinary high aspect ratios. 



3.3.1. Fiber drawing 

Special fibers with a single, centered hole with diameters down to d = 400 nm are drawn 
by the method of PCF drawing. These fibers are used as matrix for nanowire arrays. The 
drawing technique shall be briefly explained^. 

The PCFs that are used for the samples are made using a stacking and drawing tech- 
nique. The first stage is to produce a preform. This macroscopic array of capillaries with 
the same hole geometry that the final fiber is stacked by silica glass capillaries that are 
previously drawn from larger glass tubes to the desired diameters. These capillaries are 
stacked horizontally in a suitably shaped jig. The finished stack is inserted into a jacketing 

7 Nikon Eclipse LV100 microscope with a maximum magnification of lOOx. 

8 Zeiss Gemini Nvision 40 CrossBeam. The measurements were conducted with the kind help of Helga 
Hussy and Daniel Plofi 

9 An in depth review of the applications, properties and fabrication of photonic crystal fibers can be found 
in the reviews of Russell (781 [7!| . 
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Figure 3.6.: Schematic drawing of the fiber drawing and the hollow fiber filling process for sample 
fabrication. The left drawing shows a principle sketch of the drawing process for the 
fabrication of PCF or single hole capillary fiber. Inset B shows the high pressure filling 
process inside an induction over to fabricate the wires inside the silica fiber matrix. 



glass tube. The whole arrangement is mounted in the preform feed of a drawing tower 
(fig. [321(a)) . 

The preform feed slowly moves the preform into a resistance furnace which heats the 
silica glass up to temperatures near the softening temperature (tab. 13.1]) . Therefore the 
silica glass preform softens. Due to the drawing force from a cane pulling capstan the sizes 
of the preform decrease whilst intercapillary holes close, which is supported by an applied 
vacuum. The resulting canes once again are inserted into a jacket tube. 

The same process is applied again to this new preform, drawing it down to outer diame- 
ters of 125 |im or in this case to 200 um. The resulting hole size depends sensitively on all 
drawing parameters and can therefore be varied over a wide range from 300 nm to several 
(xm. 

The same procedure is applied for fabricating single hole capillary fibers where stacking 
dispenses and a silica glass tube is instead drawn in two steps and inserted into jacket 
tubes to achieve the very low diameters. 

The final silica glass fibers are coated with a UV hardened polymer coating to prevent 
it from mechanically harmful influences and to support their mechanical stability. 
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Figure 3.7.: Lateral optical microscope picture of a gold filled single hole capillary fiber (1) and 
scanning electron micrographs of gold array filled photonic crystal fibers (2-4). (4) 
shows a FIB polished cleaved fiber end face. The wires extend the fiber end face after 
cleaving the fiber due to the high ductility of gold. 



3.3.2. Fabrication of wires in fibers 

As recently success on fabricating wire arrays inside silica glass photonic crystal fibers was 
reported |7ol . 84 . 63 . for the sample fabrication in the course of this work the same 
technique was applied. 

It can be seen from the values in table [3TT1 that the melting point of the applied metals is 
much lower than the softening point of silica glass. Therefore it is possible to pump molten 
silver and gold into the holes of PCF or single hole capillary. 



The fabricated fiber is first stripped of its protective coating which would start burning 
when heated to the appropriate temperatures. The fiber is then inserted into a pressure 
cell made of thick silica glass. The insertion port is well isolated to remain high pressure, 
thus the hole structure in the fiber provides the only way out of the cell. The fiber's inner 
end is heated up above the melting point of Au, respectively Ag. At the same time its 
end face dips completely into the liquid metal. When applying a high pressure inside the 
pressure cell, the metal intrudes the fiber. After cooling down the whole arrangement to 
room temperature, thin metal wires have formed inside the fiber matrix. 

Thin wires with a diameter of 400 nm are achieved to produce successfully with this 
method. It is shown from electrical conductivity measurements, that the wires inside the 
fiber are continuous for gold samples. 



It is reasonable that due to the different thermal expansion coefficients of the metals and 
the surrounding silica matrix (a ex p, tab. 13. 1 1) . tiny gaps at the material interface form under 
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certain cooldown configurations. This is proved by an optical microscope examination from 
a side view on the fiber as well as from a comparison of previous measured experimental 
results, achieved by in-fiber optical excitation of SPP on the wire arrays to computed 
theoretical results [84]. 

The fabricated fibers containing wire arrays and single wires are cleaned from remains 
from the filling process that are found on the outer surface of the fiber using a ultrasonic 
cleaning device with a cascade of extranet, acetone (CH3COCH3), isopropanol (C3H8O) 
and ethanol (C2H6O). 

During the whole filling process they have to be handled with special care. For common 
fiber applications contaminations of the outer surface of the silica fibers do usually not 
play a significant role. In the case of side scattering, the light has to pass the outer surface 
without any measurable distortions. Especially contaminations that occur before or during 
the heating process are very difficult to impossible to remove after the cooldown. However, 
it is possible to reduce those to a sufficient level. Additionally recent further enhancement 
of the filling technique supports to avoid these effects. 



A laboratory detergent to remove especially organic contaminations that is particularly used for the first 
step of sample preparation for SEM and other high vaccum applications. 
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Chapter 




Experimental setup 

"In every branch of knowledge the progress is pro- 
portional to the amount of facts on which to build, 
and therefore to the facility of obtaining data." 
(James Clerk Maxwell, 1851 ) 

4.1. General setup 




Figure 4.1.: Principle scheme of the experimental setup. Each part of the experimental setup A, B, C 
and D is explained separately. The diagram on the right depicts the scattering process 
with incident angle a. The scattering maximum cone is drawn with the scattering 
maximum angle indicated as 7' = 2a and the angle of observation as 7. 

The experimental setup is designed for measuring optically the scattering amplitude of 
a sample, either small freestanding wire or embedded micro- to nano-wire in silica. The 

All abbreviations that are used during the explanation of the setup can be found in appendix IB. II with 
a reference to the page where the abbreviation is first introduced and explained. 
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scattering amplitude (fig. 14. ip is acquired by varying the incident angle a of a collimated 
laser beam at a wavelength of 632.8 nm. The adjustable angle of observation 7 is kept 
constant for each measurement. The maximum intensity of the reflected light is observed 
at twice the incident angle 7' = 2a, as derived in sec. 12.51 

For a constant wavelength A surface plasmon modes on the cylindrical surface are ex- 
citable only for discrete incident angles a m , as shown in sec. 12.3.21 For these angles an 
abrupt coupling to the SPP and therefore a corresponding decrease of the reflectivity of 
the cylinder is anticipated. Of special interest are lower order SPP modes as these are pre- 
dicted (see sec. I4.3[) to cause more distinct dips in the reflectivity. Particularly low incident 
angles are probed, as the mode number increases with larger incident angle a. The setup 
also allows measuring larger incident angles a as long as the range of the scattering angle 
7' does not exceed the size of the measurement sensor (PD2) which is situated close to the 
scattering centre. 

The measurements are controlled by means of a personal computer with a program that 
acquires amplitude values for each angle step by step. With the developed setup it is 
possible to measure the scattered amplitude for two different, perpendicular polarizations 
simultaneously. Thus the setup provides access to the polarization dependent scattering 
amplitudes for s- and p-polarized light (sec. I2.4| ). This not only makes a faster measurement 
possible but also minimizes artifacts from remaining vibrations or small movements of the 
probe. 

For the following explanations the entire setup is divided into parts A and B for the 
optical experiment, C for the electronics and measurement devices and D for the software, 
control and data acquisition system. 

In part A the laser beam is prepared for the measurement as described in sec. 14.21 (left 
inset in fig. l4.3p . For a simultaneous measurement of both polarizations the beam is divided 
into orthogonal polarizations, each of them being chopped by a different frequency and thus 
labelled by a frequency envelope on the signal. Hence the two polarization dependent sig- 
nals can be separated again after probing. It is possible to measure the scattering response 
for each polarization as long as the two polarizations are not mixed before scattering and 
independent of possible mixing of the polarizations after the scattering process. The two 
modulated beams are superimposed again and launched into a polarization maintaining 
fiber (PM-SMF) that guides the light to part B of the setup. 

Part B covers the scattering probing itself as depicted in fig. 14.41 and is described thor- 
oughly in sec. 14.31 The superimposed orthogonal polarizations are collimated and directed 
towards the sample. The scattering wire is rotated around the vertical axis perpendicular 
to the cylinder axis (fig. I4.ip in small steps by a motorized rotation stage (motorized RS) 
that is controlled by the software as explained in part D. A reference signal from the beam 
before scattering is acquired by a photo detector (PD1), which is connected to part C of 
the setup by BNC cables. The scattering amplitude is measured by a second, amplified 
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Figure 4.2.: The schematic diagram of part A of the setup shows the light source and the polarization 
dividing beam chopping as described in sec. 14.21 The resulting beam is coupled into 
polarization maintaining fiber (PM-SMF). Part B shows the main part of the scattering 
setup as described in section l4~3l with the reference photo detector (PD1), the rotated 
sample and the main photo detector (PD2). 



photo detector (PD2) at a constant but adjustable angle of observation 7. 



Part C comprises the electronic processing of the acquired signals, the signal separation 
for the different polarizations by means of three lock-in amplifiers, the analog-to-digital 
conversion and noise suppression. The setup is described in detail in section 14.41 



In part D the computer program is presented which was developed in the proprietary 
visual programming environment National Instruments Lab View. It provides control of 
the step motor, lock-in amplifiers and measurement devices, acquiring the measurement 
signals, processing them, cancelling noise. Additionally it acquires monitoring pictures via 
a digital adjustment microscope that looks on the rotated sample from top continuously 
allowing an error analysis of the acquired measurement data. 
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Figure 4.3.: Detailed schematic diagram of part A of the setup. The beams of perpendicular po- 
larizations in the two paths are chopped and thus labelled with different frequencies 
/i = /outer and /2 = /inner by one chopper before they are superimposed. The resulting 
waveform of the total amplitude is sketched as a function of time. 
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4.2. Beam optics 

4.2.1. The light source 

A Helium-Neon lasei0 with a power of 12 mW operating at A = 632.8 nm was utilized as a 
source of monochromatic, collimated light. Its beam is linearly polarized (> 500 : 1), single 
moded (TEMoo > 99%) with a 1/e 2 diameter of 0.88 mm and a divergence of 0.92 mrad. 

As first test measurements showed, it is necessary to stabilize the laser against backscat- 
tered light from the scattering process and further optical components. Therefore an optical 
Faraday isolatoiH (01) is introduced in the optical beam path after the laser. 

The optical isolator passes light into the intended direction just with a low intrinsic loss 
of 30%. In fact, a power of 8mW of the initial 12 mW is measured to be transmitted. At 
the same time the polarization of light travelling the optical path backwards is attenuated 
by 35 -40dB. 

After the application of the optical isolator, fluctuations of the laser are well suppressed 
to a level where referencing the measured results leads to a sufficient signal-to-noise-ratio. 
The laser is operated at least some hours after switch-on to ensure that it is stabilized. 

4.2.2. Polarization labelling by beam chopping 

With the task of measuring scattering effects for the two polarizations synchronously, a 
system was set up to deliver two perpendicular polarizations in precisely one beam path, 
each of them chopped with a different frequency. For that purpose a polarization dividing 
setup in the geometry of a Mach-Zehnder interferometer is employed as depicted in fig. 14.31 

The linear polarization of the laser beam is tilted by a half-wave plate (HWPl), op- 
timized for A = 632.8 nm and is split by a polarizing beam splitter cube (PBSljfi By 
adjusting the angle of HWP1, thus it is possible to change the fraction of horizontally 
polarized light passing straight through PBS1 (path 1) and vertically polarized light that 
is by alignment reflected to an angle of 90° with respect to path 1 (path 2). In order to be 
able to shut path 1 or path 2 separately without unbalancing HWP1 an adjustable optical 
attenuator (Al, A2) is passed by each path. 

Subsequently an optical chopped is situated, chopping of both beams at once with two 
different frequencies. Therefore the chopper wheel consists of two rings, an outer ring 
formed by 6 holes for path 1 and an inner ring with 5 holes for path 2. 

lr rhorlabs HRP 120 

2 OFR/Thorlabs IO-2D-633-VLP Free-Space Isolator, 633 nm, 2 mm Aperture, isolation 35-40 dB 
3 Linos G335592000, 5 mm, design wavelength 633 nm 

4 Stanford Research Systems SR540 optical chopper system: Frequency stability 250ppm/°C, (typ.) fre- 
quency drift < 2%, 100 Hz < / < 3700 Hz, phase jitter (rms) 0.2° (50 Hz to 400 Hz) respectively 0.5° 
(400 Hz to 3.7 kHz) 
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So beam 1 is chopped by /outer an d beam 2 by /i nri er = f /outer- The chopper itself 
is controlled by a stabilized controller (250ppm/°C with a typical frequency drift < 2%). 
The chopper also measures the two frequencies via slotted optical switches on the chopper 
wheel and two optical barriers and delivers them as reference signal. 

The two chopper hole rings are arranged in such a way that the duty cycles (the fraction 
of light time to T = 1//), of the inner and of the outer ring are similar, both 50%. This 
implicates that there are phases (30 % of the time) where both beams are open at once. 
Even if the polarization of both beams is well adjusted to be perpendicular, a decrease in 
the duty cycle minimizes the time that both beams pass the chopper simultaneously. 

Such a reduction of the duty cycle is possible by applying two similar chopper wheel 
blades stacked and slightly tilted with respect to each other. This precisely reduces the 
opening aperture of each chopper hole opening in the direction of rotation. The chopper 
controller frequency measurement for one of the two rings is obfuscated due to technical 
reasons of the frequency measurement. This problem is solved by a slight modification of 
the chopper wheels was applied to enable a correct frequency measurement again. 

At the output of the Mach-Zehnder geometry the two beams from both paths are super- 
imposed by a second polarizing beam splitter (PBS2) with low loss (typical PBS extinction 
ratio > 500 : 1). Again the horizontally polarized beam path (1) passes the beam splitter 
straight through, the vertically polarized beam (path 2) is reflected. For a resulting parallel 
and precise superposition of the two beam paths an adjustment by tilting PBS1 and PBS2 
as well as the mirrors M3 and M4 is required. 

The two superimposed beams pass alignment mirrors (M5, M6) that allow tilting and 
translation of the beam and a second half- wave plate (HWP2). For the later scattering 
part a flexible adjustment of the launching position of the beam is necessary. Therefore the 
beam is coupled into three metres of bow tie polarization maintaining 633 nm single mode 
fiber (PMJj (effective mode diameter of 3.2 \im and a numerical aperture of 0.14 — 0.18). 

The second half-wave plate (HWP2) provides rotation of the two perpendicular polar- 
izations to align them parallel to the two symmetry axes of the polarization maintaining 
fiber (PM-SMF) into which the beam is coupled by means of a fiber coupling staged. 

For efficient coupling, the beam is focused by an aspheric lens (LI) with a numeric aper- 
ture of 0.5 on a straight cleaved fiber end face by a typical efficiency of 40 to 50 %. Cladding 
modes are sufficiently removed by coiling the PM-fiber with a bending radius of 15 cm. 



5 Fibercore Limited HB 600 

6 Elliot Scientific Gold Series XYZ Flexure Stage with High Precision Manual Adjusters 
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4.3. Scattering 

To perform the actual scattering process, two different physical configurations of the setup 
are applied: 

1. A 'bottom up' approach, explained in sec. 14.3.2] is the most obvious way for probing 
scattering wires that are freestanding in air. It is also the most stable, as standing 
optomechanical standard laboratory equipment is designed for minimizing mechanical 
vibrational effects. 

Also silica embedded nanowire arrays can be probed in this configuration showing 
a superposition of the scattering effects on the surface of the supporting silica fiber, 
surrounding the wires, and the effects of scattering on the wires and wire arrays 
themselves. 

2. The latter measurement can be performed without the superposed effects, thus sig- 
nificantly more advanced, with the second configuration as described in section 14.3.31 
Surrounding the wire supporting fibers with index matching liquid, which matches 
the index of refraction of silica whilst still enabling the rotation of the sample. This 
makes a more sophisticated 'top down' setup necessary with the sample hanging down 
from the motor into the liquid. In this configuration the supporting fiber becomes 
invisible as the refractive index difference of its surface to the surrounding medium 
becomes zero. 

The setup in either configuration is particularly sensitive to vibrations as a tiny movement 
(of the order of wire's size) of the scattering target can cause a significant influence on the 
scattering measurement. To reduce setup inherent noise sources it was decided not to 
incorporate active vibration damping but to use other techniques. 

A method to minimize the influence of unavoidable mechanical vibrations and to increase 
the overall effective Pointing stability is using a larger beam diameter. For larger beam 
diameters the wire is situated in the centre of a Gaussian beam intensity profile with 
increased total width. With the same Gaussian slope thus a small movement from the centre 
does not change the incident intensity as much as for a small beam diameter. However, the 
beam diameter can not be extended infinitely in reality and an increase in the diameter 
always reduces the portion of scattered light, reducing also the signal amplitude that is 
detected. 

A suitable compromise of a beam waist of the collimated Gaussian laser beam of 800 \xm 
was chosen. Further noise reduction is realized by a sophisticated measurement and acqui- 
sition scheme as described in sections 14.4.31 and 14.51 

4.3.1. Optics 

Outcoupling from the polarization maintaining fiber (PM) is provided by a second fiber 
coupling stage in front of the scattering part of the setup. The stage is situated on the 
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Figure 4.4.: The schematic diagram of part B shows the main part of the scattering setup. As 
described in sec. 14.31 the light is collimatcd and directed onto the sample with the 
incident angle a which is stepwise successively varied by a motorized rotation stage. 
The angle of observation 7 is kept constant for each measurement. 

height level of the scattering setup and is equipped with a fiber rotation mount (FRM), 
by means of which the two perpendicular superimposed polarizations are adjusted to be 
vertically and horizontally oriented. This results in scattering with incident s-, and p 
polarization with respect to the sample and its vertical rotation axis (as defined accodring 
to sec. El}. 

It is necessary to couple in both polarizations precisely according to the symmetry axes 
of the fiber, otherwise the polarizations will start mixing over the length of the fiber. 
Measurements show that a stable separation of the two polarizations of I s : I p = 200 : 1 is 
achieved with a precise adjustment of HWP2. 

In order to prevent both fiber end faces from collecting dust and being effected by am- 
bient influences they are covered and continuously purged with nitrogen. 

The beam is collimated by an aspheric lens (L2). The e~ 2 diameter of the collimated 
beam was measured by beam profiling to be D ~ 800 \±ra. A small fraction of the beam 
is then reflected by a beam splitter (BS3) (reflectivity 8%) on a photo detector (PD1) 
for the measurement of a reference signal. This coated pellicle beam splitter (BS3) is 
not polarization sensitive, minimizes ghosting and internal reflection and thus virtually 
eliminates the amount of backscattered light from the scattering process potentially hitting 
the reference sensor. Additionally, the choice of a pellicle beam splitter minimizes the 
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lateral beam offset. 

As the angle of BS3 can be adjusted, is mechanically stabilized and is chosen to be near 
to perpendicular incidence, the portion of the reflected reference beam is stable and always 
of sufficient amplitude whilst the amplitude of the main beam is not reduced significantly 

The smaller the wire diameters, the smaller are the geometrical and effective cross sec- 
tions (compare to the theoretical treatment in sec. I2.4h of the scattering target, leading to 
a decrease in the total scattered intensity. 

This can in principle be compensated by focussing the light to a smaller Gaussian beam 
waist. On the other hand this modification leads inevitably to a broadening of the real 
incident angle a ± 6. The incident angle a being best peaked (5 ~ 0) for plane waves with 
a collimated but not focussed beam. 

The latter case of perfect plane waves was treated in chapter [5] theoretically. Any broad- 
ening of the incident angle spectrum obviously leads to a broadening of the expected 
resonance dips, thus a scattering geometry without a focussing objective for the scattering 
beam was finally preferred. 



4.3.2. Free-standing scattering setup 
Sample stages and beam centering 

The scattering sample is placed on a higher level than the other beam optical components: 
it is mounted on top of a number of stages allowing a stepwise rotation for the measurement 
itself as well as a centering of the sample in all spatial directions (fig. 14.4^ . 

As the beam direction, height, and the focal point are predetermined by the, once coupled 
and roughly adjusted, invariable position of stage 2, the whole rotation setup is placed 
on a xyz-translation stage (stage 3) with micrometer precision. By means of this, the 
sample is moved to the beam position. As this stage carries a light load but a rather 
high superstructure, it is necessary to balance it out by small weight loads to suppress 
mechanical vibrations. 

On top of stage 3 the motorized precision rotation stage is placed providing variation 
of the incident angle a for the measurement process^. The rotation stage offers a rotation 
range of > 2n with a smallest specified step width of 3.5 iirad, a maximum speed of 6°/s 
with a reverse backlash of 200 pxad. This motor is programmed to operate at a maximum 
speed of 0.25 °/s, the actual speed being calculated from the step width by the control 
program to avoid the backlash error at very small step widths. 

The rotation stage provides a magnet-and-Hall-sensor referencing instrument to ensure 
a reproducible starting position with very high accuracy. The system is controlled by a 



7 Physik Instrumente (PI) M-037.DG servo motorized, referenced precision rotation stage with PI mercury 
servo motor controller, USB programmable. 
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matched dc servo controller with feedback. 

On top of the rotation stage an xy-stage with micrometer screws is mounted (stage 4), 
used for centering the scattering sample to the rotation axis of the rotation stage. This 
manual centering whilst observing the incident beam and the wire from top by means of a 
digital zoom microscope camera (see section l4.3.4p . is the first step for adjusting the sample 
and selecting the center of rotation to be situated on the sample surface with micrometre 
precision. It is subsequently followed by the adjustment by means of stage 3 to the position 
of the incident beam. 

Sample mounting 

The thin wire samples as well as the wires supported by silica fiber are permanently attached 
to black P VC holders (picture inset in fig. I4.4p . The synthetic material offers the intrinsic 
advantages of high absorption and not causing any plasmonic effects. 

The holder is designed as a V-groove with a centred circular hole of a diameter that 
is large enough, to absorb the whole extension of the incident beam also for low incident 
angles a. Thus any possible reflection that is not caused by the sample is suppressed. The 
sample is hold straight and without movement which otherwise might be caused by air 
turbulences or resonant vibrations. 

At the same time the central hole is utilized for mounting the holder as it is threaded 
inside. By attaching the exchangeable sample holder to its mount with a transparent PVC 
screw it is possible to adjust the height of the sample to the incident beam in a position of 
perpendicular incidence. The beam causes a backward scattering pattern that can be well 
observed by means of a ground glass plate. The transparent screw lights up if the beam 
passes the sample when not exactly hitting the wire's centre. 

A side effect of this sample mounting is that half of the probing beam for an incident 
angle of a ~ is blocked, absorbed and diffusely backscattered by the black holder. 

4.3.3. In-liquid setup 

For the second type of the measurement setup, the sample is placed inside a box which is 
filled with index matching liquid^ matching the refractive index of silica. Hence, probing is 
possible without additional optical effects from the surface of the supporting silica fibers. 

Very fine surcface contaminations, relicts of the filling process (described in section I3T3]) . 
cannot be removed completely and still potentially cause scattering on interface of glass to 
index matching liquid. 

The box is made of acrylic glass with laser beam windows of silica that are slightly tilted 
against each other to suppress internal reflections at the silica to air boundary as well as 
any potential etalon and cavity effects. To maintain maximum and equal intensity of both, 



Cargille Labs: Immersion liquid 19571 n = 1.4587 ± 0.0005atA = 589.3 nman<225 °C 
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s- and p-polarization nearly perpendicular incidence instead of Brewster angle incidence is 
chosen. 

The sample holder setup as described for the 'bottom up' setup in the last paragraph 
is effectively turned upside down to a 'top down' setup as gravitation shall be utilized to 
maintain the index matching liquid confined in the box. Thus, a rotation and mounting of 
the sample from the bottom is excluded. The motor and stages are operated from the top, 
the sample with the already described holders hangs into the liquid box. 

4.3.4. Visual process observation 

To allow a visual imaging observation of the scattering measurement as well as a visual 
support for the required precise adjustment of the sample relative to the incident beam, a 
camera equipped zoom microscopqj system is focussed on the sample from top (right part 
of fig. S3]). 

The microscope tubus is provided with a 20x magnifying long working distance micro- 
scope objective with a working distance of 20 mm allowing operation whilst scattering 
measurements are running. 

The microscope is equipped with a high resolution CCD video camera (UEye) that is 
programmed and read via a USB interface with proprietary low level machine commands. 
From the zoom microscope's aperture the sample is illuminated sufficiently to acquire 
pictures of the wire with magnifications that are during the measurements usually adjusted 
to 40 x, depending on the wire diameter and sample type. 

The illumination white light does not interfere with the measurement laser beam as its 
frequency, the power line frequency of 50 Hz is filtered away by lock-in amplifiers. 

4.4. Measurement and referencing 
4.4.1. Optical detection 

Two photo detectors are used (fig. 14.41 and fig. 14.51) to detect the scattering amplitude (PD2) 
and to acquire a reference signal (PD1). The latter includes all amplitude fluctuations that 
are caused before the scattering process for a normalization of the scattering signal that is 
acquired by PD2. 

Acquisition of the reference signal (PD1) 

To measure the reference signal just before the scattering process a battery reverse biased 
Si-pin photodetector PDlj^j is operated at a wavelength and at a power level (typically 



9 Navitar machine vision 12x optical zoom microscope with CCD camera. 
°Thorlabs DET36A/M - High-Speed Si Detector, 350-1100 nm, 14 ns Rise Time 
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P ~ 400 |xW) where it provides a good linear response with a spectral responsitivity of 
7£(633 nm) = 0.4 A/W. 
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Figure 4.5.: Detailed schematic diagram of the electronics and measurement devices, part C of 
the setup. Analogue electronic signal connections are indicated in black, digital signal 
connections in grey. The signal waveforms of PD1 (reference signal) and PD2 (measured 
signal) arc fed into 3 lock-in amplifiers. Two of them are digital lock-in modules. The 
resulting DC waveforms are acquired digitally and connected to a control PC via GPIB 
connections. 

The collimated beam is focussed by a lens (L2). Its focusses size hits the sensor area 
(3.6 mm x 3.6 mm) centered. The reference sensor hence provides a reliable signal propor- 
tional to the chopped beam and not affected by possible vibrations. It effectively detects 
the sum of the intensities of both polarizations a step before the scattering. 

PD1 provides the photocurrent directly out of the photo diode anode. This current is 
proportional to the incident optical power I rc f = P1Z. It is converted into a proportional 
voltage V vc { = PlZR\ 03b A by an adjustable resistance i?i oa d = 50 O ... 50 kO. 

Depending on the fiber coupling efficiency and the desired amplitude of the beam, the 
voltage can by changing i?i oa d be easily adjusted not to exceed the input voltage level of 
the lock-in amplifier whilst maximizing the used voltage range to reduce electrical noise 
and to optimally utilize the available sampling width. 

In this case it is helpful that the modulation range of the reference signal is, apart from 
the intended chopping square modulation, low as it is only caused by the fluctuations in 
the beam intensity. Typically a load resistance of i?i oa d = 50 kf2 was applied to achieve a 
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the photo detector is the bandwidth [29j] /b w = (27t/?i oa< iCd) _1 and the rise time response 
t r = 0.35//t, w to the incident signal, being the photo diode capacitance. 

The electric circuit formed by the BNC cable, the photo detector, the load resistance 
and the measurement device is nothing else than an RC-circuit with a certain inertia, 
smoothing the edges of the chopped signal waveforms. An increase of -Ri oa d thus decreases 
the bandwidth which results in smoothing waveform edges, observable on an oscilloscope. 

For the applied chopping frequencies of / « 320 Hz a sufficient response is achieved with 
the applied load resistance. 

Acquisition of the scattered signal (PD2) 

To measure the scattered amplitude a Si-pin photodetector PD2 is applied. It includes 
the same reverse biased photo diode as PD1 but with a switchable gain transimpedance 
amplifier with an operational amplifier (in fig. 14.51 schematically depicted as A). 

The transimpedance amplifier has the advantage that, in connection with the photo 
diode, it shows only a very low resistance. The photo diode therefore shows a highly linear 
response to the incident optical power. Also for high gain factors a response is achieved 
that is sufficient for the applied measurement frequencies. The BNC connection of the 
photo detector is connected in parallel to an internal 50 f2 resistance to match the 50 ft 
impedance of the BNC cables. 

PD2 is usually operated with an amplification adjusted to 20 — 30 dB. The sensor then 
shows a total responsitivity depending on the sensor responsitivity 7£(633 nm) = 0.35 A/W 
and the transimpedance gain [Q\. 



For a transimpedance amplification of 20 dB a theoretical total responsitivity of T^tot ~ 
525 ± 2% is calculated, for 30 dB a total responsitivity of T^tot ~ 16 625 ± 2%. For some 
samples, especially those with very low diameter < 20 |xm the amplification is increased 
further. The noise of the amplifier then increases also but still remains bearable. 

4.4.2. Principles of lock-in amplification 

Lock-in amplifier^ are used in a broad spectrum of applications due to their ability to 
measure signals with a high dynamic range (the ratio of noise contributions to the measure- 
ment signal) and a very low signal-to-noise ratio. They offer phase sensitive measurement 
and the separation of frequencies with a very narrow bandwidth. 

11 Thorlabs PDA 100A switchable gain, amplified silicon detector, 400 — 1100 nm 

12 Only the most important features of lock- in amplification are mentioned here. For an in depth discussion 
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Figure 4.6.: Schematic diagram of the signal processing of the lock-in amplifier. Two phase sensitive 
detectors (PSD) are displayed for measurement of amplitude (X) and quadrature (Y) 
as realized in dual channel lock-in amplifiers. The right inset displays schematically 
the processing of an edge smoothed rectangular input signal as it is typical for optical 
chopping measurements. 

In the setup for this work the lock-in amplifier's ability to separate signal contribu- 
tions with different frequency modulation and the cancelling abilities for noise of other 
frequencies are particularly utilized. The second advantage cancels also contributions from 
disturbing ambient influences caused by e.g. other experiments in the same laboratory. 
Particularly all influences from laboratory illumination or power line interferences are can- 
celled as the power line frequency of 50 Hz is filtered away very effectively. 

The operating principle of lock-in amplifiers with special attention on the features that 
are practically utilized in the setup is explained in the following (fig. 14.61 

The voltage controlled oscillator (VCO) is a sinusoidal (RS830) or rectangular (lock- 
in modules 1 and 2) waveform generator that is triggered by the frequency / s p of the 
square shape reference signal F re f (frequency locking). The signal Vvco that is generated 
by this internal oscillator can be adjusted in phase to fit the measurement signal V- m which 
is important for the next step (phase locking). The RS830 can automatically determine the 
phase mismatch, thus it is also applied (preliminarily triggered with the other frequency 
f p ) for the determination of the phase for p polarization that is then manually configured 
into lock-in module 2. 

In this experiment no substantial phase shift over the measurement time is expected, 
hence for any phase shift far enough from n/2 a constant signal would be expected. Thus 
a correct phase locking to 0° simply provides a maximum output signal which is useful to 
minimize later electronic noise contributions. 
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The phase sensitive detectors (PSD, also called demodulators or mixers) are the 
core parts of any lock-in- amplifier. Their function is in principle the multiplication of the 
reference oscillator signal Vvco with the measurement signal V\ n that already is roughly 
cleaned from other frequency components by a cascade of adjustable notch filters. 

A short mathematical explanation may clear the functional principle of the PSD. As- 
suming a noise-free sinusoidal signal voltage with angular frequency u at the input V- m = 
Acos(tot). The reference signal V rc f with the same frequency and a respective phase shift 
of —ip triggers the VCO to generate a signal with zero phase shift Vyco = B cos(uit + ip). 

The output of PSDx is then obtained by a multiplication of both signals: 



For a constant amplitude B (VCO) and a constant frequency uj, thus the output is a 
DC signal, that is proportional to A (signal), the cosine of the phase difference ip and 
modulated with 2ut, which physically means that the first harmonic component remains. 

Technically the PSD consists of signal splitting, inversion and triggered gating elements. 
It is followed by a low pass filter providing time averaging for an adjustable time scale, called 
integration time (which is usually chosen to be 300 ms) and whose purpose is to eliminate 
the 2u>t modulation. For a lock-in amplifier this averaging is thus intrinsically necessary, 
revealing the major disadvantage of the instrument: Lock-in amplifiers are comparably 
slow, becoming slower for higher measurement accuracy. 

In all applied lock-in amplifiers the PSDs are realized as digital electronic circuits with 
the main advantage of real linear multiplication with a high dynamic reserve and good 
harmonic rejection. 

The dynamic reserve of a lock-in amplifier is the adjustable ratio of the largest tol- 
erable noise signal to the full scale signal, quantitatively expressed in dB. At the lock-in 
amplifier the dynamic reserve at which it may be operated without an input overload can 
be selected. Higher dynamic reserve allows measuring of signals with a large noise con- 
tribution at the expense of a low gain and a/d conversion noise. Thus when applying a 
lock-in amplifier always a compromise has to be found. 

The signal outputs of the lock-in amplifiers deliver dc signals, proportional to the 
incoming ac signal amplitude at the correct frequency and phase. X out represents the 
amplitude component with 0° phase shift. This is passed for the reference signal at the 
output of the lock-in module 1 only as this module has a single PSD design. RS830 and 
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lock-in module 2 deliver also the quadrature amplitude Y ou t (90° phase shift) and the radius 
R of the vector in a X-Y- plane. The latter value, R is acquired for both scattering signals 
for further data processing as this has the advantage to be intrinsically independent of any 
phase shifts. 

4.4.3. Polarization dependent lock-in detection and noise handling 

The two polarizations (s and p) being scattered by the sample are detected by only one 
photo detector (PD2). This 'one photo detector' principle makes it necessary to separate 
the two signals after detection. The separation is realized using of lock-in amplifiers that 
are locked to the chopping 'frequency labels' of the two different polarizations (sec. l4.2."2"|) . 

Electronic separation of the measurement signals 

The application of 'one photo detector to detect both polarizations' has the major ad- 
vantage that both polarizations can be observed in exactly one angle of observation syn- 
chronously. A little spatial offset would already cause a large offset in the angle of observa- 
tion 7 as the sensor PD2 is placed very near to the sample to capture a maximum scattered 
light. 

On the other hand the conversion from optical power to electric current by the photo 
diode at this point already has reduced one information dimension, only detecting the 
intensity I = \E a + E p \ 2 and not the two dimensional (in the plane of the sensor) electric 
(or magnetic) field vectors E StP . 

The clearer the polarizations of the two beams, chopped with / s and / p are separated 
(i.e. E s _L Ep), the less interference between both polarizations occurs. In the case of a 
non zero parallel component of the beam portions from path 1 and path 2, interference is 
unavoidably caused by slight optical path length differences, already at the order of the 
wavelength of light, in the Mach-Zehnder like chopping setup. 

These interference effects would affect the measurement accuracy thus reducing the ac- 
quired information of the whole measurement irreversibly. As a separation of the two 
polarizations of 200 : 1 is practically measured to be achieved with the chopping setup for 
precise PM-fiber coupling (as mentioned in section H~3]) the interference effects between the 
two polarizations are low enough to allow the desired measurements with high accuracy. 
This interference reduction is supported by a reduction of the chopping duty cycle (as 
described in sec. 14.2.2ft 

In principle, optical effects that mixed the polarization of both 'beams' at the scattering 
process whilst not destroying the coherence of the light could lead to the same interference 
noise. This is not expected to happen from theoretical assumptions. It is shown by oscillo- 
scope measurement of the signal waveform (fig. l4.7p that with a precise adjustment of the 
beams and optical components interference effects in fact do not play a significant role. 
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Figure 4.7.: By means of a digital oscilloscope acquired typical signal waveform (PD2) during a 
measurement (compare to figure I4.3[) . The two amplitudes of the s-polarized (green 
line, 0) and the p-polarized (blue line, <-*•) beam portion are indicated. The noise level 
on each single amplitude arc of the same order as the noise amplitude with both beam 
portions at once, indicating that interbeam-intcrfcrcncc is very low. The remaining 
noise is mainly due to laser noise and is mostly cancelled by reference division. 

Scattering signal PD2 S 

PD2 is connected to the single voltage signal input of a dual channel digital DSP lock-in- 
amplifier SR830 El (fig. l4.5p . The lock-in amplifier's internal sinus oscillator is externally 
synchronized to the (reference) frequency of the s-polarized beam (/inner = fs) which is 
delivered from the frequency generator and chopper control. 

The two superposed signals of the s- and the p-polarized beam portions are of the same 
amplitude. Additional amplitude noise caused by interference, detection or optomechanical 
parts of the setup on the PD2 signal (fig. 14.71) cause a total signal-to-noise-ratio (SNR) of 
the order of 1 x 10 3 — 1 x 10 4 . The laser noise contribution to the total noise is higher but 
still much lower than the signal amplitude. 

Therefore the lock-in amplifier can be operated with low dynamic reserve (sec. I4.4."2|) . 
thus allowing a higher measurement accuracy as the a/d noise and the dc gain noise de- 
er easd^ 4 !. The low dynamic reserve setting of the SR830 changes for different gain settings. 

13 Stanford Research Systems SR830, for detail it is referred to appendix [Cl 

14 The SR830's low noise adjustment still allows a relatively high dynamic range even in low dynamic range 
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Not discussing this dependency in detail it shall be mentioned that the sensitivities (maxi- 
mum full scale signal detection amplitude) of the SR830 are usually, depending on a good 
combination of the PD2 gain factor and the SR830 and lock-in module 2 gain factors, 
chosen to be 200 mV, 500 mV, 1 V. With an adjustment of 'low dynamic reserve' this quan- 
titatively results in a dynamic reserve of 4 dB, 6 dB, dB. 

The SR830 is programmed and the measurement data is acquired via its GPIB PC 
interface. 

Scattering signal PD2 P 

The photo detector PD2 for the scattering signal is not only connected to the SR830 but 
also to the signal input of the lock-in module 2l£| (fig !4.5p . This instrument is an integrated 
lock-in module that in principle operates equal to the SR830 with two (amplitude and 
quadrature) phase sensitive detectors (PSDs). 

The main difference is that it lacks the ability to automatically lock the phase which has 
to be adjusted manually. The vector radius output R is obtained as for the RS830. Thus 
referencing the phase might seem unnecessary which is in fact not true as for a wrong phase 
adjustment the simpler (compared to the SR830) square signal VCO (not a sine generator) 
could cause a/d signal mismatch artifacts. 

The lock-in module 2 is on its reference signal input connected to the frequency generator 
and chopper control's output for the /outer frequency that labels the p-polarization beam 
portion. 

As the lock-in module lacks an integrated a/d converter it is connected to a programmable 
digital precision multimeteillfl which is programmed via its GPIB interface. The full range 
R-output amplitude of the lock-in module 1 is maximized to 10 V. The detection precision 
of the multimeter is adjusted respectively to achieve a maximum a/d conversion precision, 
minimum a/d discrimination artifacts and to minimize electronic noise. 

Reference signal PD1 S 

The output of the photo detector PD1 is connected to the signal input of the lock-in 
amplifier module 1 (fig. 23]). This instrument is predominantly identical to the lock-in 
amplifier module 2 but has only one PSD and provides only the phase 0° amplitude X out 
as an output, not the quadrature Y out or the radius R ou t- 

The reference signal input of the lock-in amplifier module 1 is connected to the fre- 
quency generator and chopper control frequency output /i nnC r together with the SR830, for 
frequency locking on the s-polarized beam portion. The lock-in module 1 in the beginning 

mode of operation, compared to analogue lock-in amplifiers as it is constructed to minimize the total 
noise level and a further decrease in the dynamic range would not lead to a significantly further decreased 
noise level due to the digital detection process. 

15 Femto LIA-MVD-200L - Digital dual channel lock- in module. 

16 Agilent 34410A, 6.5 digits digital multimeter 
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of the measurement cycle is manually adjusted to the phase of the frequency reference for 
s-polarization. This phase does in fact not change over the measurement and not signifi- 
cantly over different measurements, also for this signal always the maximum 0° amplitude 
is achieved to acquire with the same advantages as mentioned before. 

It was chosen only to acquire the reference signal for one of the two chopping frequencies 
thus beam polarizations (s) as preliminary measurements showed that for an optimized 
coupling the noise of the beam signals of both polarizations does not differ significantly. 

The SR830 is in addition to its main purpose also utilized as an a/d converter with the 
ability to intrinsically acquire samples exactly synchronized to the SR830 lock-in sampling, 
when correctly programmed. Therefore it supplies external DC signal input channels which 
though provide only a a/d conversion resolution of 16 bit. 

After some tests it was decided not to use these input channels for the measurement of 
the scattering signal PD2 p as the electronic noise contribution distorted the measurement 
signals too severely. Nevertheless, for the scattering signal a very high dynamic range over 
the whole measurement, from signal to maximum range amplitude has to be detected, 
which is not necessary for the reference signal PD1 S . 

The reference signal is expected to contain only the noise fluctuations that it is intended 
to provide for referencing, no additional amplitude drift over the measurement as it is 
modulated on the scattered signals. Therefore the a/d convertors of the SR830 with the 
reference signal are always driven at maximum input amplitude, there providing the full 
16 bit resolution for the whole measurement process which is sufficient for the reference 
signal. 

The SR830 provides the a/d convertor signal for PD1 S on the second of two measurement 
channels via the GPIB programming interface. 

4.5. Data acquisition and control software 

The great advantage of lock-in amplifier detection is its low measurement noise level and 
its capability to separate electric signals of different frequency (sec. I4.4."3]) . On the other 
hand the operating principle of a lock-in amplifier makes measurements comparably slow. 
Additionally a requirement for a large a and amplitude resolution requires a high number 
of measurement steps. A complete measurement over an incident angle range a = 0° . . . 10° 
with an angle resolution of 0.01° takes more than 3 to 5 hours. 

At least these reasons motivate a complete automation of the experimental measurement 
including a computer assisted processing of the acquired measurement values and several 
error corrections. The software to control and automate the experimental measurements is 
written in the visual programming language National Instruments Lab View 8.2. For precise 
synchronization the communication with the measurement instruments was programmed as 
low level command functions, embedded in the visual programming language. A simplified 
schematic of the software elements is depicted in fig. 14.81 
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Figure 4.8.: Schematic diagram of the measurement and control software, part D of the experimental 
setup. 



The interface of the software part D to the electronic measurement instruments of the 
setup C (see section l4.4.3p is realized by the digital data communication protocol GPIB. Via 
this protocol the RS830 digital lock-in amplifier and the digital multimeter are controlled 
and their data is read and transferred for processing into the PC and the control software. 
The motorized precision rotation stage (setup part B, fig. l4.4p is controlled via a USB 
connection. It is programmed with a proprietary low level command language as well as 
the control camera of the zoom microscope that observes the scattering process from top 
(see also fig. I4.4p . 

The whole software is modularized and includes precise error handling and diagnostics. 
The measurement results are saved in reduced, human readable format on hard disk. 



4.5. Data acquisition and control software 
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Initialization and referencing 

To start the program, the measurement parameters are input by the user. It is necessary 
to drive the motorized rotation stage to its reproducible reference start point, to partly 
manually adjust the sample in the beam and to drive the stage to its starting position 17 !. 
After a change of the sample slight differences in the holders can change the a = 0° 
position. This has to be adjusted by observation of the scattering pattern and the measured 
amplitude as well as the controlled zoom microscope camera. 

The measurement parameters are the lock-in sensitivity, the sampling rate for the syn- 
chronous acquisition of the measurement data, the angle a range of measurement and the 
number of steps to measure, respectively the angle a resolution. 

Finally all instruments are set to the measurement parameters and prepared for acquisi- 
tion. All further communication with the instruments is software triggered. 

Measurement loop 

The measurement loop itself consists of three subsequent steps. First the sample is driven to 
the desired incident angle a. In the next step synchronized scattering amplitude detection 
takes place. In the third step the observation microscope camera is triggered to acquire a 
control image. 

The motor for each position reports back the exact real position that will be used for 
further steps and for the data acquisition including backlash correction. The position is 
measured and calculated by the servo motor's parameters and counter. 

For the acquisition of the scattering data it is necessary to achieve a precise synchroniza- 
tion of the different measurement devices (SR830 and digital multimeter). These commands 
are implemented on a basic hardware command level with correction factors that are the 
result of test measurements. 

With a lock-in hardware integration time of usually 300 ms and a sampling rate of usu- 
ally 128 Hz, adjusted to fit to the measured low frequency noise contributions, up to 10 s 
data for the scattered s-polarization via the SR830 (PD2 S ) is acquired. For the scattered 
p-polarization the multimeter is read synchronously to achieve a PD2 P waveform. The 
reference signal PD1 S waveform (as described in section I4.4.3H is acquired via the internal 
a/d converter unit of the SR830. 

Signal processing 

All three waveforms show the noise modulation of the scattered beam. The main portion 
of the signal noise is cancelled or at least reduced significantly by the lock-in amplifiers. 
Remaining noise on the lock-in output signals is dominated by low frequency contributors 
that modulate the square chopped signal amplitude and pass the frequency filters. The 



The rotation stage is equipped with a Hall sensor for referencing and counts its motor position from that 
point of reference. 
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4. Experimental setup 



so called laser residual intensity noise (RIN) [29( is the dominant contributor to the total 
noise level of the scattering beam and is far above the shot noise level. HeNe lasers usually 
show a very low phase noise but this does neither way influence the measurement in this 
experimental setup. 

As the total noise modulation is not assumed to simply disappear by a time integration it 
is an increase in real measurement precision to divide the signals of s(t) := PD2 s (t) /PDl s (t) 
and p(t) := PD2 p (£)/PDl s (i) for the waveforms pointwise and after that calculating the 
time averaged mean values s and p. 

By displaying the acquired waveforms from the different photo detectors the coincidence 
of the low frequency (1 — 100 Hz) common noise can clearly be identified. Before the refer- 
ence division signal-to-noise-ratios of 1 x 10 2 to 1 x 10 3 are observed, after the pointwise 
reference division this increases to 1 x 10 3 to 2 x 10 5 always depending on the conditions, 
amplitudes and gain factors of the current measurement. 

To simplify a later data analysis also the mean values of the original scattered amplitudes 
and the mean value of the reference signal are processed. To achieve a measure of the 
waveform fluctuations over time, the standard deviations a of all values are calculated and 
bundled simultaneously. As the light source fluctuations are divided out, the standard 
deviation of the divided values should be lower than the original ones (usually <7 SjP ~ 
1 x 10" 4 ). 

After the measurement of the scattering amplitudes a reference picture is automatically 
taken from top view by the microscope to allow later manual crosschecking of observed 
peaks to possible not intended reasons. 

The calculated scattering amplitudes (s, p) are displayed over the incident angle a and 
stored. The loop continues its functions with increasing a until the measurement is finished. 

Closing the measurement 

After finishing the measurement loop the systems are driven back to their initial configura- 
tion. All acquired data along with possible error messages and the measurement parameters 
are saved in a measurement file. 
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Theoretical results 



5.1. Numerical computation of the dispersion model 

The model for SPP modes on cylindrical samples, derived in section 12.3.21 is used to com- 
puta3 the dispersion of SPP (eq. I2.35|) and the angles of incidence (a) for exciting modes 
(m) (eq. 12.37]) . The computations are accomplished for the experimentally relevant material 
configurations of silver (Ag) or gold (Au) in silica (Si02 glass) or air. 

For the dispersion of the dielectric permittivity of silica glass, experimental values (see 
tab. 13.11 fig. I3.ip and the Sellmeier equation are applied. For the optical properties ofgold, 
experimental data is used (fig. I3.2p . Silver is modelled by the Drude model (fig. I3.3|l 2 l 

In figure [5TT] the dispersion of SPP modes on the surface of a silver cylinder with diameter 
20 |xm is shown, calculated from equation (|2 .35[) . For a better visibility only modes m = 
1, . . . , 101 in steps of 10 are shown. 

The plot can be directly compared to the case of planar SPP (fig l2.3p . It shall be pointed 
out that, in comparison to the undamped dispersion (fjg. l2.3[) . also for planar SPP a finite 
maximum wave number (3 is approached at the surface plasmon frequency wsp if rea l 
materials are considered. 

In this case, including damping, the quasibound, leaky part of the dispersion relation 
in the region wsp < <jJ < lo p is not forbidden any more. This "backbending effect" can be 
observed in all dispersion curves which consider losse^l. The effect can also be observed in 
figures 15.11 and 15.21 

The depicted curves correspond well to the theoretical dispersion relations that were 



calculated by Pfeiffer et al. 67j, as long as an undamped Drude model is assumed instead 



1 These computations are implemented in computer programs, written in the proprietary programming 

language Matlab. 
2 For details about the material properties see sec. 13.11 

3 To avoid redundancy, for the case of planar SPP it is here not depicted explicitely. It can be easily found 
in common textbooks, e.g. Maier and Atwater [H2], p. 29. 
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Figure 5.1.: Dispersion of SPP on a silver cylinder of D = 20 |a.m, surrounded by air. Displayed are 
the dispersions of modes m = 1 (in black) to to = 101 (in grey) in steps of 10. The 
left figure shows the dispersion in terms of the normalized SPP wave vector (3cq/uj v 
(cq. I2.35p . The inset is a magnification for the wavelength range around the laser 
wavelength (A = 632.8 nm). The right figure depicts the incident angles for excitation 
of the same modes (eq. I2.37|) . 



5.1. Numerical computation of the dispersion model 
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wavelength A 
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Figure 5.2.: Dispersion of SPP on a gold cylinder ffig. !3.2[) of D — 20 urn, embedded in silica 
(fig. EHJ • Displayed arc the dispersions of modes m = 1 (in black) to m = 151 (in 
grey) in steps of 10. In both figures the wavelength of the applied laser A = 632.8 nm 
is indicated by a dotted line. The left figure shows the dispersion in terms of the 
normalized SPP wave vector /3co/w p (eq. I2.35p . The right figure depicts the incident 
angles for excitation of the same modes (eq. I2.37[) . 
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angle of incidence a in ° angle of incidence a in 



Figure 5.3.: Dispersion of SPP modes on cylinders of gold in silica (left) and silver in air (right) 
over the radius and incident angle. The excitation angles are calculated (eq. I2.35|) for 
a wavelength of A = 632.8 nm. Displayed are the modes m = 1, 2, 12, 22, . . . 

of the complex dielectric permittivity. Pfeiffer et al. presented an exact solution of radiative 
SPP modes below lo p on the surfaces of metal cylinders for this simplified case. 

The inset (fig. 15.11 left) displays the dispersion in the wavelength region of experimental 
interest. The laser wavelength at A = 632.8 nm is indicated as a horizontal dashed line. For 
the experimentally relevant wavelength regime it can be clearly seen that the dispersion 
curves of the modes lie above the light line. Thus an excitation of these radiative SPP 
modes is possible. 

For a gold wire in silica, as depicted in figure [53| the dispersion is slightly more compli- 
cated than for a silver wire due to the interband absorption in the visible regime (compare 
fig. I3.2p . It is qualitatively the same as for a freestanding gold wire, only the remaining 
edge from the forbidden region in the undamped planar case smoothes out even further, 
which also smoothes the gap in the plot over the angles of excitation for further modes of 
lower order. Therefore a second plot for the air-gold configuration is not presented here. 

In the limit of very large diameters the dispersion of the cylindrical SPP approaches the 
dispersion of planar SPP as presented in section I2.2( for an infinite diameter, a cylinder is 
a single planar interface. This behaviour is obviously confirmed analytically, referring to 
eq. (I2T35]) . 



5.2. Comparison to an exact analytic solution 



S3 



For both configurations depicted in this section, the angles for excitation of the modes in 
a wide range of radii (0.1 \xm < R < 20 (Jin) are calculated, applying eq. (|2.37|1 . The range 
of radii as well as the wavelength (A = 632.8 nm) are chosen according to the experiment. 
The resulting angles a are depicted in figure 15.31 

It is apparent from this chart that the SPP modes for gold cylinders in silica spread 
further than those for silver cylinders in air. The modes of low order can be seen to 
converge for increasing radius at low a in all cases. Higher order modes converge with 
decreasing radius towards a = 90°. From the model it can be therefore expected that for 
low a, thus in the case of striking incidence, for sufficiently large radius a number of distinct 
low order SPP excitations is observed whose spacing is very small, followed for larger a by 
an increasing mode spacing. This prediction is done under the assumption that the model 
is valid in the region of interest. 



5.2. Comparison to an exact analytic solution 



Schmidt and Russell [83j have reported a different approach to determine the dispersion 
relation of long-range spiralling SPP-modes on metallic nanowires. Applying the theory 
usually used to calculate guided modes in optical wave guides j9l[ to the structure of 
metallic wires in dielectrics, an exact solution for the dispersion of different kinds of modes 
on the surface is obtained. 

In detail, the approach starts from solving Maxwell's equations in cylindrical coordinates, 
applying boundary conditions on the metal-to-dielectric interfaces. The resulting dispersion 
relation can be written as 18311: 



(qfyi + <7fVo)0i<?oV'i + £o<?fVo) - m 2 n 2 m (e - e ; ) 2 = 

where m is the azimuthal mode order. The radial wavevectors outside q Q and inside the 
cylinder q\ and the abbreviations ip Q and ip\ are defined as follows: 

K m+1 (k q o R) J m+1 (k qiR) 

tfj = m - k q o R — ipi = m - koqiR—— — 

K m (k q o R) K m (k qiR) 

q Q = ±\A4i - £ o <?i = V £ i - n m 

J m (x) being the Bessel function of first kind, K m (x) being the modified Bessel function 
of kind (both of order m) . n m represents the effective refractive index of mode m as also 
described in the model (eq. (I2.36|) ). In the mentioned publication the validity of the model 
whose results are discussed in sec. 15.11 is cross checked with computed dispersion results 
from this exact approach. It is found to agree increasingly well for higher order modes. 

For a determination of the range of validity of the model in the relevant radius and 
wavelength ranges for radiative modes (above the silica respectively vacuum light line), the 
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Figure 5.4.: Comparison of the dispersion of cylindrical surface modes from the derived model (lines) 
to the dispersion from an exact solution of the guided modes (dashed). Calculated arc 
modes m = 5 ... 12 on a gold cylinder (R = 2 ... 10 urn) in air. 



solutions for the appropriate modes are searched and the corresponding dispersions are 
calculated0. 

The obtained effective refractive indices n m are then converted to excitation angles ana- 
logue to eq. (|2.37p . 

A comparison of the resulting angles a for radii 1 urn < R < 3.5 |J.m from the model of 
spiralling SPP and the exact solution is depicted in figure [5741 

The model shows good agreement with the exact solutions for lower radii per mode, 
respectively for higher a for all modes. The deviation of both solutions increases with lower 
incident angles. Particularly for small angles from the graphs in fig. 15.41 it is apparent that 
the dispersions do not coincide. 



5.3. Numerical computation of the scattering coefficients 



For this analysis again it is important how the material properties are described. Martinos 
and Economou [54| point out that modelling the excitation of surface plasmons in scat- 



4 These calculations were kindly contributed by Dr. Markus Schmidt. 



5.3. Numerical computation of the scattering coefficients 
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Figure 5.5.: Absorption efficiency for different angles and diameters on a single of gold and silver 
in air. The left graphs (1,3,5) show the amplitude of the scattering efficiency Q a b s as 
a colourplot over a variation of incident angle (a) and radius (i?, in (i.m), calculated 
from eq. for a wavelength of A = 632.8 nm. The graphs in the right column (2,4,6) 
overlay the angles of SPP excitation as calculated from the model (eq. 12.371 fig. I5.3[) 
over each of the left graphs for comparison. Only <5 a b s for p-polarization is shown, as 
the s-polarization does not exhibit a resonance structure. 
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tering theory makes it necessary to explicitly introduce damping of the electromagnetic 
waves inside the material. This is especially the case as many theoretical treatments of 
the problem omit losses and additionally assume a dielectric permittivity of the metal of 
^metai = ~~ 1 resulting in strongest possible SPP resonance, although it is not experimentally 
applicable. 

The application of only the real part of e(u), thus omitting losses, in fact simplifies the 
numerical solution of the scattering problem significantly: Complex values as arguments of 
the Bessel functions, Hankel functions and their derivatives in the equation system (|2.4ip 
result in magnitudes for physically real parameters lead to problems in the numerical 
representation of 32 bit computer programs^. This was an even more significant challenge 
in earlier studies, as described by Shah j8^|. 

In the described calculations the wavelength dependent intrinsic dissipative properties 
of the material are entirely included by applying complex dielectric permittivities e(to). 

The first step to solve this problem in this work was to derive a linear equation system 
and the corresponding solution matrix (|2.4ip . This matrix is numerically solved instead of 
computing the coefficients a ns , a np , b ns , b n ^{ from equations (12.44h . which is a much more 
efficient approach and provides a higher precision of computation. 

The wire radius R enters the arguments of the Bessel functions r/,£ fj2.42l2.43p . respec- 
tively of the whole expansion. It is shown that the higher R is, the more orders n of the 
expansion are necessary, to include all nonzero coefficients. The expansion is shown to 
converge for orders m > 15 for small radii {R < 3 (xm) respectively small angles . For a 
computation of all incident angles 0° < a < 90°, (e.g. in fig.s I5.5|5.7p . the order m of 
the expansion has to be increased to m < 220. A routine is implemented in the program 
to determine the necessary number of orders and to indicate results that cause erroneous 
results. 

The range of radii for certain angles can due to this method be enlarged, compared to 
a computation from equations (|2.44p . It is therefore possible to calculate the scattering 
efficiencies for radii up R = 20 [im e.g. for silver and gold wires with surrounding air. 

Applying this method, extinction and scattering efficiencies Q ex t and Q sca are calculated 
from eq. (|2.47p . Eventually the absorption efficiency is obtained from eq. 12.461 

The curves of the absorption efficiencies Q a bs ( n g- I5.7|) in the case of s-polarization have 
substantially different shapes than those for p-polarization: 

For s-polarization Qabs increases with a (fig 15. 7p . For the wavelength that is also experi- 
mentally applied (middle gray, dashed, 632.8 nm), the curve increases monotonic and 
no resonance dips are observed. For shorter frequencies, this behaviour stays the 
same. 

5 An overview over the Bessel functions for complex arguments is compiled in appendix lA.il 

6 To achieve the scattering efficiencies, obviously only the coefficients a, b for the fields outside the cylinder 

are necessary. / and g for the expansion of the fields inside the cylinder are only used for the computation 

of the field distribution as shown in section [57^1 



5.3. Numerical computation of the scattering coefficients 
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For p-polarization the slope of Q a bs remains positive. Different to s-polarization, a modu- 
lation of the absorption efficiency over a is clearly indicated. 

For longer wavelengths (black, 1200 nm) the absorption efficiency shows anomalous be- 
haviour. For both polarizations the curves are modulated by an additional fluctuation. 
Particularly in s-polarization this results in a modulation significantly visible in the region 
of small incident angles. It is not assumed that this superimposed effect is caused by the 
interband transitions of gold, as (see fig. I3.2p these lie at much shorter wavelengths. A 
precise examination of the effect could reveal interesting results but is postponed as the 
focus of this work is the examination at A = 632.8 nm. 

The left plot in figure 15.51 displays the absorption efficiency Q a bs for a silver wire in air 

(So = ifl 

Comparable plots of the computed scattering cross section over the wavelength A and a 
variation of the refractive index for coated cylinders with perpendicular incidence with the 
objective to examine finite wavelength cloaking have been recently published Nicorovici 



et al. [60] . 



Martinos and Economou [54| demonstrated that the imaginary part of the eigenfrequen- 
cies of the surface plasmon modes increases for higher mode order n, which causes the 
sharpness of the peaks in the absorption coefficients Q a bs to decrease with n. This theoret- 
ical prediction can also be clearly seen in the graphs for p-polarization in figure For 
increasing a, thus higher mode orders, the dips smooth out. 



5.3.1. Inherence of SPP resonances 

It can be shown that the denominators of equations (I2.44b( I2.44b|) coincide morphologically 
with the resonance angles a calculated from the exact solution of the dispersion relation 
of SPP modes on the cylinder (as presented in sec. I5.2|) . 

This observation suggests clearly that the derived scattering model accords to the descrip- 
tion of SPP modes on the cylinder surface. The same is therefore true for the approximated 
model in its region of coincidence with the exact solution. 

Englman and Ruppin already discussed a theoretical comparison of an exact solution 
for the dispersion of modes on ionic crystals of cylindrical shape to the scattering patterns 
calculated analogue to Q&bs- Although the size of their particles lies clearly in the regime 
much smaller than A, in the limiting case their result also shows a coincidence of the 
resonances in the scattering theory and those from an exact solution of the mode dispersion. 



Recently She et al. 89( have reported a theoretical comparison of exact classical scat- 



tering theory to analytic formulae for plasmonic resonances on small wires {D <C A) with 



7 It shall be particularly mentioned that the blue (= 0) area for radii R > 17 urn at small angles is a 
numerical artefact. The coefficients for these radii cannot be obtained from the applied theory, for the 
aforementioned numerical reasons and are therefore set to zero. 
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low dissipation. Their results also indicate an interconnection as the one proposed here, 
although their calculations were again only computed for ideal e\ = — 1 + i8 with very 
small dissipation, which is then in detail transferred to very thin, coated Ag cylinders. The 
detailed results therefore cannot be compared. However, the resulting resonance behaviour 
can be compared to the one achieved in this work and shows coincidence. 

As the comparison of an exact solution for modes of low order (sec. 15.2ft indicated a 
rather good agreement, at least for sufficiently large R and a, the calculated R over a for 
wires of gold (2,4) and silver (6) are overlayed to the amplitude of Qabs (right plots (2,4,6) 
in fig. l5.5p . Particularly for large angles of incidence a a matching of both theories can 
be observed. However, for incident angles a < 20°, the calculated modes do not coincide 
with the absorption efficiencies. As the range of a in the experiments is particularly chosen 
in this range, therefore it is decided to compare the measured results primarily with the 
computed results for Qabs- 



5.4. Numerical computation of the scattered field distribution 

Although the scattering coefficients (sec. I5.3|) allow a good insight into the angle depen- 
dency of the total scattering properties of the cylinder, a more detailed analysis of the 
scattered field distributions is worth being done. The field distributions and the energy 
flux, represented by the Poynting vector offer near field information about the scattering 
processes in contrast to the far field scattering properties of the scattering efficiency (as 
introduced in section |2.4.6|) . 

Scattered field distributions, the Poynting vector field distribution and Poynting vector 
field bifurcations for the simplified case of normal incidence (a = 7t/2) on small (D <C A) 
coated cylinders were recently studied in detail by Luk'yanchuk and Ternovsky [491 ] . 

The field distribution (E(r)) rtVtZ is calculated for all three cylindrical coordinates by 
solving matrix (I2.4ip analogue to the algorithm for the computation of the absorption 
coefficients Q a bs- F° r the field distribution inside the cylinder E\, not only the scattering 
coefficients a n ,b n but also f n ,9n are determined. 

After a transformation from cylindrical coordinates to carthesian coordinates, the field 
distribution is plotted for an area chosen to show the wire surface with high resolution. 

Figure [531 depicts the absolute field distributions (\E\) exemplarily for incident plane 
waves with a = 45° (top six pictures) and for normal incidence (a = 90°, bottom six 
pictures) as well as the projections of the Poynting vector. Each row of pictures shows the 
three perpendicular spatial planes in cartesian coordinates with respect to the cylinder and 
the plane of rotation, as illustrated in fig. 12.51 

In all pictures, the plane waves approach the cylinder from the right side. 

Examining the field distributions over a variation of the incident angle reveals the for- 
mation of spots of field enhancement on the wire surface. The insets in figure 15.71 show 
that for both polarizations such field enhancement effects occur. It is assured by several 
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Figure 5.6.: Field distribution and energy flux, calculated for a = 45° and a = 90°. Plane waves 
with A = 632.8 nm strike a single gold wire D — 1 u.m in air. The absolute value of the 
electric field vector (|-E|) and the values and directions of the Poynting vector (S) arc 
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checks of order of the field expansion and from checks of the continuity of the tangential 
field components that the observed effects are not computing artefacts. 

An striking effect is that also in the case of s-polarization, for increasing a, on the non 
illuminated side of the wire for a particular choice of parameters, spots occur. These can 
be seen in figure 15.71 whilst for the parameters in figure 15.61 the effect cannot be observed. 
It seems that this effect for the s-polarization is caused by localized surface plasmon (LSP) 
excitation, which cannot yet be proved due to restricted computation time. 

Particularly for p-polarization, for which from the model SPP excitation is expected, 
not only with increasing a discrete field spots on the surface occur. From a fine resolved 
simulation of field plots over a, a periodic fluctuation of those is observed, which can be 
correlated to the modulation of the absorption efficiencies as depicted in figure 15.71 

The Poynting vectors for p-polarization indicate an energy flow along the surface of 
the wire. Different to the simplified schematic in figure I2.4[ the beam, in the theoretical 
calculation the plane wave, extends not only over one point on the top or bottom half of 
a horizontal cylinder but hits its full lateral extension. This means that if the resonance 
condition is fulfilled, SPP modes into both, counter-propagating directions are expected to 
be excited. 

An outstanding feature of a comparison of the field distributions is the formation of a 
"shadow" behind the wire in the case of s-polarization, but the emission of radiation from 
the wire in p-polarization that causes no "shadow" to arise. 



Comparing the depicted resonance behaviour to recent publications [89l . |49j , analogies 
in the field distribution are observed. For these treatments of small cylinders with a very 
thin coating, as in the case described here, with higher mode order that is achieved by 
a variation of the frequency of incident light, discrete modes arise, which are represented 
by additional nodes of field accumulations on the first and second surface (outer coating 
interface and inner coating interface). 

The pattern of field distributions, observed in this work is a bit more complicated, al- 
though the complification of a surrounding thin layer is not treated. Additional modes 
are observed to occur at increasing a and a reorganization of the existing modes is seen 
in the final field plots. New modes arise at the upper and lower edge of the wire in a 
perpendicular cut as depicted in the right field plots and the spacing of the existing mode 
spots decreases. 

Due to the much higher range if h mc that is accessible by a variation of a compared with 
the reported, common variation of u in the mentioned publications, a much higher mode 
order is reached already for relatively low incident angles. 



5.4. Numerical computation of the scattered field distribution 



91 



4.:, 



Qabs 



4 - 



3.5 - 



2.5 



2 - 



4.5 - 



0.5 



-3 -1-1 



-3 -1 _ 1 3 




s-polarization 400 lira 
p-polarization 400 ran 
s-polarization 632.8 ran 
p-polarization 632.8 ran 
s-polarization 1200 ran 
p-polarization 1200 ran 




40 



20 



30 



40 



50 



incident angle! ol in ° 



DID 



I -El 



i -3 -1 j. 1 3 

a 3 = 18.1° p-pol. 



Figure 5.7.: Absorption efficiency Q a b s of a gold wire in silica (i? = 3|xm). Plot for s- and p- 
polarization at the frequencies A = 400 nm (light gray), A = 632.8 nrn (middle gray) 
and A = 1200nm (black) and resonance field (\E\, arbitrary units) and Poynting vec- 
tor distributions for 632. 8nm at angles of 7° (1), 9° (2), 18° (3) for p- (top) and 
s-polarization (bottom). Spatial units in urn. 
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Experimental results and discussion 



"Nur Bcharrung fiihrt zum Zicl, 
Nur die Fiillc fiihrt zur Klarhcit, 
und im Abgrund liegt die Wahrheit." 
(Friedrich Schiller) 

Using the experimental setup as described in ch. 0] series of measurements of scattering 
amplitudes are accomplished on wires (ch.[3|) of gold and silver. These two materials are 
chosen due to their good processing properties (sec. 13. 1 [) and as for these two materials the 
filling method for fibers in PCF and single hole capillary (sec. I3.3| ) shows the best results. 

An additional major advantage of wires of gold, fabricated by the sputtering technique is 
its extraordinary durability against surface degradation. In fact, a comparison of the results 
confirms that the measurement of SPP resonances on silver wires becomes increasingly 
difficult much faster than for gold wires. 

A sufficient number of free-standing sputtered samples shows distinct, reproducible dips 
for the theoretically predicted angles of SPP excitation. Therefore a concise comparison 
of the acquired data with theory is presented in sec. 16.21 Good accordance between both 
is indicated. Therefore the first proof of radiative SPP on wires for nonperpendicular inci- 
dence is reported to be accomplished successfully. 

Also samples of small silver and gold wires and wire arrays in specially drawn silica fiber 
(sec. I3.3|) were fabricated successfully that are by electric conductivity and optical lateral 
examination proved to be continuous. The measurements of these samples with both 
methods (sees 14^2] I4.3.3P turned out to be practically much more difficult. Scattering on 
single wires was first tried. The effective cross section of these samples was measured to be 
too small to obtain results with dips in the scattered intensity. Possible effects from SPP 
resonances could not be separated from measurement noise. 

Therefore the scattering properties of samples of wire arrays in PCF were measured. The 
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6. Experimental results and discussion 





Figure 6.1.: Scattering pattern of a single, free-standing silver wire with a diameter of D = 20 urn. 

The displayed plane is the plane of the sensor (PD2), in which a ground glass screen 
is positioned for these pictures. The patterns are acquired under a variation of the 
incident angle a by photography. 



complicated resulting patterns could not yet be correlated to theoretical effects. In this 
case coupling between the modes on adjacent wires occurs, causing a modification of the 
resonance properties. A further optimization of the setup to allow a successfull measure- 
ment and the derivation of a corresponding theory (which could base on the propositions 
of Pitarke et al. [68l j) were not accomplished in the course of this thesis. 



6.1. Observed scattering pattern 

In accordance to the theoretical scattering pattern (sec. I2.5|) . the experimental setup ex- 
hibits circular rings in the experimentally accessed range of incident angles (0° < a < 20°). 
This observation is in accordance to the low eccentricity of the theoretically predicted 
ellipses due to low a. 

Displayed (fig. l6.ip are the observed half rings, as the samples are mounted on a holder 
(see sec. 14.3.2)) that intentionally blocks the incident laser beam. This causes also blocking 
of the left half of the scattering patterns. 

The observed intensity distribution suggests that the main portion of light reflected by 
the cylinder, is radiated into a relatively small total solid angle and that it is strongly 
peaked on the surface of the scattered ellipses. This supports a measurement of scattering 
with a collimated beam, which is the best real approximation of plane waves: the relatively 
small portion of the beam that is scattered provides only a low intensity. This, on the 
other hand can be exploited well. Only small intensity is lost by scattering into the whole 
sphere surrounding the sample. 



6.2. Measurement of resonance angles 
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Figure 6.2.: Example for a typical intensity profile for scattering on a gold wire with R = 20.3 u\m 
in s- and p-polarization. This measurement is labelled as (f) in the summary and in 

fig. ES 

For a variation of a with the detector (PD2) (sec. \AA\ at a fixed angle of observation 7 
(fig. I4.ip therefore, the maximum intensity step by step passes the aperture of the sensor. 

Best results are achieved with a choice of 10° < 7 < 20° for measurements with incident 
angles 0° < a < 10°. As PD2 is placed very near to the sample and offers a comparably large 
aperture, followed by the focussing optics, it is possible to catch the scattered maximum 
already for low a and track it for the mentioned angle range with high sensitivity. 

6.2. Measurement of resonance angles 
Measurement of the scattered intensity distribution 

For the reported measurements continuously improved techniques are applied to produce 
the sputtered gold wires. A significant enhancement is observed by an optimization of the 
sputtering process (sec. l3.2p . 

For each of the samples a number of measurements are carried out. The measurement 
parameters are chosen for a range of 0° < a < 10°. The strongly peaked intensity maximum 
that can be clearly seen in figure 16.11 for each measurement travels through the aperture 
of the sensor, leading to a maximum at 7 = 2a (see fig I4.ip . This intensity profile is 
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6. Experimental results and discussion 




Figure 6.3.: Graph of the scattering results for gold wires of different radii in air. The horizontal 
lines indicate the radius for each measured wire. The symbols indicate the mean of the 
rcproducibly measured dip positions. The data is ovcrlaycd over simulated absorption 
efficiencies from fig. 15. 51 

modulated by the absorption profile of the sample, which is to observe. 

All measurements are carried out with an angle resolution of at least Aa < 0.01° with an 
optimized resolution of the acquired signal amplitudes. An increase of the angle resolution 
always causes an increase in the total time necessary for a complete measurement. 

Depending on the measurement conditions, a typical automated measurement takes two 
to four hours. Smaller diameters of the wire samples reduce the intensity of scattered 
light as the effective scattering cross section decreases. The consequently reduced signal- 
to-noise-level requires longer integration times (sec. I4.4."3|) . For very small wire diameters 
nevertheless the noise on the resulting measurement curves cannot be avoided to increase. 

For a reliable identification of the small dips caused by SPP excitation, at least a ratio 
of the SPP dip-strength to noise of 5 is found to be required. 

The referenced signals for both incident polarizations (s and p) acquired simultaneously 
(sec. I4.4|) are compared for each measurement. 

An exemplary curve of the reflected intensity over a variation of the incident angle a is 
depicted in figure [6\2l The angle of observation for this measurement is chosen as 7 = 16°, 
thus at an incident angle of a = 8° the maximum intensity (see fig. 16. If) passes the centre 
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Figure 6.4.: Graph of the scattering results for silver wires of different radii in air. The horizontal 
lines indicate the radius for each measured wire. The symbols indicate the mean of the 
rcproducibly measured dip positions. The data is overlayed over simulated absorption 
efficiencies from fig. 15. 51 

of the sensor aperture. 

The referenced intensity of s- and p-polarization in normalized, arbitrary units is dis- 
played together with the ratio of p- over s- polarization. The ratio reveals the first dip 
(indicated with an arrow) to be measured only in p-polarization. The dips at higher a are 
measured in both polarizations but reproducible for repeated measurements. 



Analysis and reproducibility 

As it is shown in fig. 15.51 a theoretical comparison unveils that the SPP model is not a 
good approximation particularly for low incident angles. Therefore for figures (I6.3)6.4[) it 
is decided to compare the experimental data directly to the results of the scattering theory. 
This is legitimate for a proof of SPP as it is shown in sec. 15.3.11 that the scattering model 
inherently comprises the resonance condition for SPP modes. 

Each measurement of a sample showing resonance dips is repeated at least twice to repro- 
duce the observed resonance angles. After at least three coinciding repeated measurements 
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6. Experimental results and discussion 



sample 


R in pm 


identified dips 


std. dev. \<j\ (in °) 


a 


6.5 


3 


0.12 


b 


9.7 


3 


0.05 


c 


14.3 


3 


0.08 


d 


15.2 


4 


0.15 


c 


15.3 


3 


0.21 


f 


20.3 


3 


0.20 



Tabic 6.1.: Measured results for samples of sputtered gold wire of different radii. For each sample 
the measurement was repeated at least twice. The measured number of dips and the 
mean standard deviation for all dips over the measurements are indicated and plotted 
in fig. 16.31 



sample 


R in pm 


identified dips 


std. dev. <7 (in °) 


a 


10.5 


3 


0.23 


b 


11.9 


3 


0.15 


c 


12.6 


3 


0.20 


d 


16.15 


3 


0.09 



Table 6.2.: Measured results for samples of sputtered silver wire of different radii. For each sample 
the measurement was repeated at least twice. The measured number of dips and the 
mean standard deviation for all dips over the measurements arc indicated and plotted 
in fig. 16.31 

of at least three dips the sample is included in the overview of results for gold (tab. 16.1]) 
respectively for silver (tab. I6.2|) . For a comparison of the resonance dips to the calculated 
absorption efficiencies of gold (fig. l6.3p . respectively silver (fig. I6.4j) the mean value of the 
angles of these reproduced dips is calculated and the standard deviation is indicated. 

In figures (16. 3( \QAh it can be seen that not for every probe subsequent dips are indicated. 
This is due to a strict evaluation if a resonance dip is at least reproduced for three times. 

Both graphs indicate that at least no resonances were measured in the maxima of Q a b s , 
but even the dip spacing corresponds well to the calculated distribution. 

Limiting conditions 

In the experimental data not only dips at the predicted angles are observed but also a fine 
superstructure occurs. 

However, it is not possible to exclude an additional source for absorption experimentally 
that might explain those: Coupling caused by surface roughness. In experiments to probe 
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planar SPP excitations the systematic creation of such surface roughness in the size of 
10 nm to 100 nm is applied to obtain coupling conditions. 

The surfaces of the probes, applied in this work are intended to be as smooth as possible. 
By inspectation with atomic force microscopy (AFM) during the optimization process for 
sample fabrication (sec. [3]) a surface roughness of 5 nm for the desired thickness of the layer 
was verified for the probes applied. 

This technique should minimize effects from coupling due to surface roughness but cannot 
cancel them. It is difficult to predict the resulting structure of excitations by random surface 
structuring but it can be assumed that a portion of the observed fluctuations, at least of 
those that are observed repeatedly by a repetition for one probe, is caused by this effect. 

Another particularly observed effect is degradation of measurement features during the 
measurements. It is not expected due to the exceptionally low coupling efficiency to the 
SPP modes that this effect is caused by damage from energy accumulation in resonance as 
it happens for in-fiber excitation experiments on nanowires. 

A much simpler reason is revealed by a comparison of samples of gold and silver. For 
the measurements of silver wires it is observed that the distinction of the resonance dips 
decreases more rapidly with time after fabrication, instead the number and density of noise 
dips increases. This effect is attributed to the formation of an oxide layer on the silver 
surface. 

This does not happen for gold samples due to its significantly lower chemical reactivity. 
Actually, for gold samples measurements over a longer period still exhibit the observed 
features much better. 

The only possible way to prevent this effect for free-standing samples of metal is to encase 
the setup by a vacuum chamber. Another possibility, extending applying the proposed 
fabrication technique could be to use the new possibility of sputtering deposition of silica 
glass. A thin layer of silica, surrounding the sputtered gold wire would prevent the metal 
layer from any degradation effects caused by oxidization. Of course, as a side effect it 
would significantly change the properties of the SPP, but not in a way that would not be 
explainable by the theory that is presented in this work. 

However, embedded metal nanowires in silica fiber intrinsically do not suffer this limitation. 

Additional to the difference, observed between samples of gold and silver, also for gold 
a degradation of the resonance dips over time is observed. As a reason for that the un- 
avoidable adsorption of dust can be identified which is in most cases not removable without 
destroying the sample. 

Fluctuations between curves obtained for a repeated measurement for one sample can 
be caused by the limitations of the measurement itself as well as by the assembly of small 
particles on the surface of the wires over time. 
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6. Experimental results and discussion 



Comparison to theory 

In figures 16.31 and 16.41 for both materials good accordance of the measured resonances to 
the theoretical absorption efficiency Qabs is clearly indicated. 

From theory (sec. [2]) it is anticipated that only for one polarization (p, TM, which means 
perpendicular to the plane of rotation of the wire) SPP excitation occurs. Nevertheless the 
measured data reveals dips in the scattered intensity in both polarizations for most of the 
samples. For no sample an undisturbed curve for s-polarization is obtained. 

It was experimentally checked that this effect is no measurement artefact as the separa- 
tion of both simultaneously measured signals is of the magnitude of 1 : 200. Mixing and 
interference between the two polarization channels can therefore not be the reason. The 
resonance dips physically occur for the indicated polarizations, contradictory to theory. 



A comparison to data, obtained by Miziumski [571] . [581] . [591] who showed a proof of concept 
of a scattering setup for the excitation of surface plasmons on fairly large (0 ~ 50 iim)0 
tapered glass cylinders with a smooth evaporated surface coating of aluminium and silver 
of unknown but sufficient thickness, reveals that his measurement for low angles matches 
the resonance angles from the derived model for the planar SPP approximation. 

But in fact, the comparison of this model to the exact solution and to the theoretical 
absorption efficiencies, as shown in sec. 15.31 do explicitly not coincide with the model for 
the range of low incident angles that is also probed by Miziumski. 



*A measurement for smaller wire diameters was proposed, which were not feasible at that time. In those 
experiments no polarization dependency was measured. 



101 



7 



Chapter 

Conclusion and outlook 



This thesis shows a coherent treatise of the excitation of surface plasmon polaritons on 
small wires in theory and in experiment. The measured experimental results predominantly 
agree with the modelled results from theory. The results of this work are summarized by 
formulating answers to the main questions that are established in chapter [TJ 

"What is the range of validity of the proposed model for the propa- 
gation of radiative SPPs on metal wires?" 



The model of spiralling surface plasmon polaritons (sec. 12 .3.2|) 84l . l58l . l59l . 1671 ] bases on 
the assumption that plasmons on a cylinder can be approximated by planar SPP (sec. 12.21) . 
The approximation appears intuitively to apply better for large wire diameters. This 
assumption is proved and compared to an exact solution of Maxwell's equations for SPPs 
on cylinders (sec. I5.2|) . Both models show good agreement for sufficiently large incident 
angles (a > 5°). For small incident angles it is not valid as can be clearly seen (fig. I5.4p . 

Quantitatively the model matches the exact solution better with increasing radii (R), 
higher order modes (m) and larger angle of incidence. 

"Does a description of the scattering of light on metal wires comprise 
SPP resonances?" 

For comparison, a theoretical description of the scattering process was derived from 
classical Mie theory. The inherence of SPP resonances in this theory is proved (sec. I5.3.T]) . 
Comparing the angles of SPP excitation to calculated absorption efficiencies Q a bs from this 
theory (sec. I5.4|) . good agreement for large angles of incidence but decreasing agreement in 
the range of a < 20° is observed (sec. I5.3.T]) . 

"Does an experimental scattering examination of exciting SPPs on 
small metal wires reveal resonances that match the predicted dis- 
persion relation in a range of wave vectors that is not feasible by 
in-fiber coupling?" 
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7. Conclusion and outlook 



An experimental setup was developed for the polarization dependent examination of 
scattering (sec.H]) on samples of small wires of gold and silver (sec. E]). 

A low signal-to-noise-ratio was achieved by applying sophisticated noise reduction tech- 
niques (sec. I4.4."3|) allowing the observation even of small SPP resonance effects. 

Measurements of the scattering amplitude over a variation of the incident angle exhibit 
reproducible dips at those angles, predicted from Mie scattering theory and SPP dispersion 
theory. A comparison of the experimental results for scattering on single wires of silver 
and gold (6.5 [im < R < 20.3 |J.m) shows good agreement to the calculated absorption 
efficiencies (sec. 16. 2|) . A polarization dependency of the observed resonance dips could 
not be consistently observed. In most measurements the resonances that were identified 
occurred in both polarizations. 

"Is it possible to relate the observations from scattering to a near 
field description?" 

For an in depth understanding of the occurring processes of SPP excitation by scat- 
tering, the theoretical approach from Mie theory for far field scattering efficiencies was 
extended (sec. I2.4."4"j) . With the obtained theory the near field distribution for s- and p-po- 
larization was examined for the occurence of resonance features (sec. 15.41) . 

The prediction from the SPP model, that analogue to planar SPPs, resonances only 
occur for p-polarization (TM) is proved by effects in the near field pattern and in the far 
field absorption efficiency. Arising field enhancements in p-polarization are observed in the 
vicinity of resonance angles in the near field distribution, which are identified as SPPs. The 
amplitude of these fluctuates smoothly with resonant versus not resonant incident angles. 

Additionally field enhancements in s-polarization on the surface of the wire can be seen, 
being probably identified as nonretarded localized surface plasmons (LSP). In contrast the 
amplitude of SPP, no oscillating fluctuation in the strength of these field enhancements 
can be observed. 

In the near field distribution for both different effects (s-polarization, LSPs, p-polariza- 
tion, SPPs) a remarkable difference was observed which accords to the general expectation 
for LSPs and radiative SPPs. In the vicinity of some |J.m behind the irradiated wire, in 
s-polarization a "shadow" is observed, in p-polarization this is not seen, which matches the 
interpretation as radiative SPP intrinsically radiate their energy away into the surrounding 
dielectric medium. 

A further examination of the excitation of SPPs on small metal wires still promises inter- 
esting insights. Especially as the theoretical treatment in this area of physics currently is 
ahead its experimental examination and most theoretical studies do not match the experi- 
mental requirements of materials and technical feasibility. Additionally some observations 
during the work on this thesis could not yet be explored in depth and some developed 
techniques promise further, exciting applications and further developments: 
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The calculated near field patterns could be observed by extending the developed experi- 
mental setup by a scanning near field microscope (SNOM). A comparable approach for a 
planar glass-gold-air system was recently reported (33|. Particularly the detailed exami- 
nation of the observed shadowing effect and differences in s- and p-polarization that could 
not be measured by the far field measurement technique, applied in this work, could reveal 
interesting conclusions. 

Not presented in this work are the obtained computational results for perpendicular 
incidence on a single nanowire with a variation of the wavelength instead of the incident 
angle. These computational results promise interesting experimental results. An experi- 
ment would be the complementary to the experiments carried out in this thesis. A possible 
realization of the setup could utilize a supercontinuum light source and a monochromator 
for scanning the frequency and measuring the absorption. The simulated results propose 
an interesting characteristic not only for the reflected but also for the transmitted light 
from a collimated beam. An experimental observation of the scattering effects on large 
(0 > 100 pin) dielectric cylinders for perpendicular light incidence and varied angle of 
observation and wavelength was reported by Abushagur and George [lj]. A comparable 
experiment on large diameter (0 ~ 120 (J.m) copper and brass cylinders in the infrared 
range was performed by Cohen et al. These results could be easily achieved and by 
far extended by the proposed setup. 

The techniques that were applied in this work offer the possibility of the fabrication of 
very small wires as well as dielectric wires with a very thin coating for the experimental 
testing of recently published theoretical studies. This makes it possible to probe cylindrical 
glass-metal-air configurations if a further optimization of the sputtering parameters could 
be achieved. 

Measuring experimentally the scattering properties of arrays of gold nanowires during 
the works for this thesis revealed more complicated resonance behaviour. This has to be 
theoretically and experimentally examined in detail. Particularly it is expected that in 
this case coupling between the wires occurs, changing the total dispersion relation. This 
cannot simply be modelled by classical scattering theory any more. An approach to solve 
this problem could follow the propositions of Pitarke et al. [68| . 

The fabrication technique of the samples in this works could be extended to create wires 
with a nonuniform diameter over the length. Coupling light into the metal surface at a point 
of smaller diameter, followed by a cavity of large diameter and subsequently a decrease in 
diameter could result in a cavity resonance and allow an insight into the absorption of 
travelling SPP. A modification of the surface by application of focussed ion beam (FIB) 
or other microstructuring techniques also promises an interesting modification of the SPP 
guiding properties of the wires. 

During the works on this thesis a proof of concept for a technique of side polishing the 
fiber embedded wire samples was demonstrated. This approach could allow coupling to 
these very thin metal wires without disturbance due to the bent surface of the surrounding 
fiber. Additionally ATR (attenuated total reflection) coupling using prism coupling setups 
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7. Conclusion and outlook 



could provide access to the guided SPP modes on the wires that were not possible to probe 
by the techniques applied in this thesis. 
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Appendix 



A 



Mathematical and electromagnetic 
foundations 



A.l. The Bessel differential equation and its solutions 

Due to the special importance of Bessel functions for optics in general and for the deriva- 
tions of section 15.31 in specific, some important properties of the Bessel functions shall be 
assembled here. A more in depth information about the Bessel functions and differential 
equation can be found in Stegun and Abramowitz 94], Bateman and Erdelyi 0. 



The Bessel differential equation is the solution of the scalar wave equation 12.231 in cylin- 
drical coordinates (eq. 12.4.3]) . It can be expressed as: 



d 2 _ d 



d^ z + ^ z + (l2 -" 2)z = 



Where Z are Bessel functions. The Bessel functions of first Zn\x) = J n { x ) and sec- 

(2) 

ond kind Z\ (x) = Y n (x) (the latter also called Neumann function) are the two linearly 
independent solutions to the Bessel differential equation. For integral orders n these two 
solutions are depicted in figure IATT1 over real arguments. 

The first derivatives of the Bessel functions of type v, Zn^ (x) are given by the relatiorQ: 

Z^' = \{Z { :l l (x)-Z^ +1 {x)) 



See Stegun and Abramowitz [gj for further details. 
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A. Mathematical and electromagnetic foundations 




Figure A.I.: Plot of the Bessel functions of first J n (x) and second kind Y n (x) for real arguments x 
and different orders n = . . . 5. 

A third and a fourth solution of the Bessel differential equation can be constructed from 
J n (x) and Y n (x), which are called Bessel functions of third and fourth kind or Hankel 
functions of first Hn\x) and second kind Hn\x): 

HP (x) = J n {x)+ iY n (x) (x) = J n (x)-iY n (x) 

For the calculation of the coefficients as it is accomplished in sections 15.31 and 15.41 for 
the expansion of the fields in vector cylindrical harmonics (equations 12.391 and 12.4.3)) it is 
important to estimate the number of orders of Bessel functions that is necessary to achieve 
correct results. 

Analytically a summation of inf orders is possible, in numerics the number of orders is 
crucial as the values of the Bessel functions have to be calculated for each order thus the 
time for computation scales linearly with the number of orders. 

For typical arguments of the Bessel function H n (x) of x = i] = 0.248 + 0.6477i (for a 
gold wire of diameter R = 20 \xm with surrounding silica at an incident angle of a = 45°) 
a plot over the orders n = —4 ... 4 is depicted in fig. IA.2I 



A. 2. The Poynting vector of plane waves 
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Figure A. 2.: Plot of the Hankel function of first kind and order H\(x) for real arguments and a 
plot over integral and fractional orders n = —4- • -4 for a typical argument H n (x — 
£ = 0.248 + 0.6477i) as it occurs for a R = 20 fim diameter Au wire in silica with an 
incident angle of 45°. 

A. 2. The Poynting vector of plane waves 

The intensity mathematically is derived from the nonobservable Poynting vector S of an 
electromagnetic wave. 



I{r) = \S{r) 



S(r) = 1/2 Re (E(r) x H{r)*) 



The Poynting vector includes additional information about the direction of energy flux. 
In general S is a complex vector, whose real and imaginary parts represent the physical 
energy flux respectively an analogon to the idle power. In this thesis the Poynting vector 
is defined as only the real part. 

For a plane wave (see section |2.1.4[ p. ITS]) , whose propagation direction is || e the Poynt- 
ing vector S can be written in terms of only the electric or the magnetic field. Therefore 
it is possible [47J, p. 7 ff] to describe the electric field respectively the magnetic field as: 



E 



e e 



ex H 



H = ±J^exE 
Mo A* 



(A.l) 



Thus E, H and e form an orthogonal trihedron, describing a wave which propagates in 
a homogeneous medium. The Poynting vector of a plane wave can then be written as: 
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A. Mathematical and electromagnetic foundations 



S = E x H* 



e e / V V MoA* 



= - ((e x ff) • E*) e + ((e x if) • e) £?* 
The second term is zero. The first term remains and can be converted with equation IA.1I 



to: 



S = ±J^\E\ 2 e 
Mo A* 



S = ±Jf^\Hfe 

V £o£ 

Equation IA.2I is also true for H (equation IA.2|) H with equation IA.ll Additionally for a 
plane wave we achieve the basic correlation: 

fiofi\H\ 2 = e^e\E\ 2 

The intensity of an electromagnetic wave (electromagnetic power per area), incident on 
a surface perpendicular to S is defined as: 



2 As (a x b) ■ c = — (a x c) ■ b 
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Appendix 



B 



Conventions 



B.l. Mathematical and physical conventions 

Throughout this thesis attention was paid to consistent mathematical, physical and exper- 
imental conventions. These are listed in the following lines. 



Complex relative dielectric permittivity 

In this work e = e r + isi is defined 
dimensionless. This nomenclature 
was consistently chosen for reasons of 
mathematical convenience. It shall be 
pointed out that this is not the con- 
vention of some theoretical textbooks, 
including Born and Wolf and Jack- 



son 



31] 



The conversion between both systems 
is simply e = e a bs/eo- 



Also Bohren and Huffman [12j in the 
beginning of his derivations applies 
their absolute dielectric permittivity 
nomenclature, whereas he later omits 
this completely as he expresses all re- 
lations in terms of k which is indepen- 
dent of the nomenclature of e. There- 
fore the later derivations of the scatter- 
ing theory in Bohren and Huffman [l^] 
are nevertheless directly comparable. 



Complex quantities: The real and imagi- 
nary of a complex quantity are writ- 
ten in the most widespread notation 
as x = Re(x) + ilm(x). 

Constants , physical as well as mathemati- 
cal are not typeset italic but upright, 
including n. 

Equalities are distinguished into nondirec- 
tional =, logically directional defini- 
tions := and =:, defining the quantity 
on the : marked sinde. Approximate 
equality is indicated with w, combina- 
tions are = is only used 
to emphasize a relation in contrast to 
another relation. 

Labelling of equations is limited to those, 
which are referred to in later text or 
equations. 

Vectors are typeset boldface italic a, latin 
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letters as well as greek letters, follow- 
ing the angloamerican standard print 
notation. Unit vectors are typeset e. 

Vector products are typeset a ■ (3. 

Vector differential operators are noted as 



V which is a vector whose entries are 
the first derivatives into each spatial di- 
mension. The Laplace differential op- 
erator, A A = V • (V • A) which is built 
on that convention, shall not be con- 
fused with (V • V)A. 



B.2. Nomenclature 

The applied nomenclature in the theoretical derivations and the abbreviations in the exper- 
imental sections are chosen consistently. The following overview lists all those, mentioning 
the pages in parentheses where they are first mentioned. The list starts with mathematical 
and physical conventions and constants, that are mainly used in formulae and texts, that 
describe derivations. It continues with a list of the abbreviations, mainly used for the 
explanation of the experimental setup. 



Nomenclature and constants 

a Incident angle of light with re- 
spect to cylinder respectively wire 
axis. (p. 28) 

(3 Propagation constant of surface plas- 
mons. (3 m represents discrete prop- 
agation constants, /3 represents in 
this geometry the propagation vec- 
tor, (p. 27) 

Co The vacuum speed of light is a con- 
stant, defined by the special the- 
ory of relativity, c = ^//j.qEo -1 = 
2.997924 58 x 10 8 m/s. (p. 8) 

e Complex relative dielectric permittiv- 
ity e = E\ + \£2, which is dimension- 
less. Expressed in terms of the refrac- 
tive index e\ = n 2 — k 2 and £2 = 2nn 
(See also EE]), (p. 8) 

£0 The dielectric permittivity 
(also called 'diel. con- 
stant') of vacuum is eo = 
8.8542 x 10" 12 (As)/(Vm) (p. 8) 



E[ nc The incident electric field vec- 
tor, (p. 33) 

E Q ,E\ Electric field outside (o) respectively 
inside (i) the cylinder, (p. 33) 

E sca The scattered electric field vec- 
tor, (p. 33) 

H Q ,Hi Magnetic field outside (o) respec- 
tively inside (i) the cylinder, (p. 33) 

7 Angle of observation, angle of the 
measurement photo detector with re- 
spect to the incident beam. (p. 68) 

H[ nc The incident magnetic field vec- 
tor, (p. 33) 

H sca _ The scattered magnetic field vec- 
tor, (p. 33) 

Hn\x),Hn\x) Hankel function of type 1, 
respectively type 2 and order n of the 
argument x. (p. 100) 

J n (x) Bessel function of type 1 and order n 
of the argument x. (p. 99) 



B.2. Nomenclature 



111 



k(oj) Propagation vector of an electro- 
magnetic wave. The absolute 
value is called propagation constant 
k(u). (p. 18) 

k Q , ki Wave vector of an electromagnetic 
wave in the medium outside the 
cylinder (o) or inside the cylinder 
(i). (P- 35) 

A Wavelength in m, usually indicated 
as [A] = 1 nm. The relation be- 
tween wavelength and frequency is 
A = c/u. (p. 28) 

H Complex relative magnetic perme- 
ability, in dimensionless nomencla- 
ture. Wavelength dependent mate- 
rial constant. Defined consistent to 
eas/i = /iabs/W (p- 8) 

/xo Complex magnetic permeability u. = 
4nx 10~ 7 (Vs)/(Am). (p. 8) 

M Macroscopic magnetization of a 
medium, (p. 9) 

n Complex index of refraction, defined 
by n = n r + in n r = yfejl. (p. 10) 

v Frequency [v] = s -1 = Hz. (p. 13) 

0(x n ) Landau symbol. A function / E 
0(x n ) has a dependency on variable 
x of order n. Following the definition 
precisely the Landau symbol O de- 
fines an asymptotic upper bound for 
the function /. (p. 17) 

uj The angular frequency u = 2nv [oj] = 
s- 1 = Hz. (p. 13) 

uj p Plasma frequency of a mate- 
rial, (p. 13) 



P Macroscopic electric polarization of a 
medium, (p. 9) 

r Spacial coordinate vector, usually 
three dimensional. May be defined 
different in different coordinate sys- 
tems, (x, y, z) in cartesian coordi- 
nates, (r, if, z) in spherical cylindrical 
coordinates, (p. 32) 

R\oad Adjustable load resistance for 
PD1. (p. 78) 

a The specific conductivity of a mate- 
rial, (p. 9) 

v Velocity, particularly the phase veloc- 
ity of an electromagnetic wave, v = 
~ k - (P- 18) 

Y n (x) Bessel function of type 2 (Neumann 
function) and order n of the argu- 
ment x. (p. 99) 

Zn^ (x) Bessel function of type v and order 
n of the argument x. (p. 99) 

Abbreviations and setup components 

Al,A2 Adjustable optical attenua- 
tors (p. 71) 

BNC Bayonet Neill Concelman - coax- 
ial connectors and cables with an 
impedance of 50 $7. (p. 71) 

BS3 Coated pellicle beam splitter (p. 74) 

/ Technical frequency, particularly the 
signal envelope frequencies, unit 
[/] = a -i = Hz. (p. 71) 

FRM Fiber rotation mount (p. 73) 
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GPIB General Purpose Interface Bus/IEEE- 
488. Digital communication interface 
system, applied to control the mea- 
surement and motor systems, (p. 83) 

HWP1,HWP2 Half wave plates with a de- 
sign wavelength of 633 nm (p. 71) 

L1,L2 Aspheric biconvex lenses (p. 72) 

L2 Focussing lens (p. 78) 

Laser HeNe laser (p. 69) 

M1-M6 Mirrors (p. 72) 

microscope CCD microscope with ad- 
justable zoom tubus, (p. 76) 

motorized RS Motorized precision rotation 
stage, programmable, (p. 75) 

01 Optical isolator (p. 70) 

PBS1,PBS2 Polarizing beam split- 
ters (p. 71) 

PD2 p Experimentally acquired signal por- 
tion for the p-polarization from photo 
detector PD2 (p. 83) 

PD1 S Experimentally acquired signal por- 
tion for the s-polarization from photo 
detector PD1, the reference sig- 
nal, (p. 84) 



PD2 S Experimentally acquired signal por- 
tion for the s-polarization from photo 
detector PD2 (p. 83) 

PM Polarization maintaining single mode 
fiber, (p. 72) 

PD1 Photo detector 1 (p. 77) 

PD2 Photo detector 2 (p. 78) 

PSD Phase sensitive detector. Part of the 
lock-in amplifiers, (p. 80) 

SPP Surface plasmon polariton, sur- 
face bound electromagnetic exci- 
tation (retarded) of the electron 
plasma, (p. 21) 

SP Surface plasmon, surface bound 
electromagnetic excitation (non- 
retarded) of the electron plasma, (p. 21) 

SR540 Optical chopper system (p. 71) 

stagel,stage2 Fiber coupling stage with mi- 
crometer screws, (p. 72) 

stage3 xyz-translation stage with mi- 
crometer screws and long traverse 
path. (p. 74) 

stage4 xy translation stage, (p. 75) 

VCO Voltage controlled oscillator. Part of 
the lock-in amplifiers (p. 80) 
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Applied instrumentation 



Sample manufacturing and preparation 

Tapering rig Explained schematically in section 13.2.11 For tapering silica glass fibers to 
the desired diameters. 

Sputtering system Emtech K575X turbo single magnetron sputter coater, peltier coold. 
For coating tapered fibers with thin layers of metals. 

Sputtering system AJA, ATC orion series uhv, water cooled multiple magnetron and RF 
sputter coater. For coating tapered fibers with thin layers of metals. 

Optical microscope Nicon Eclipse LV100 optical bright field/dark field digital microscope 
with a maximum magnification of 100 x. 

Scanning electron microscope Hitachi S-4800 field emission sem. 

Scanning electron microscope / focused ion beam system Zeiss field emission SEM/FIB 
Gemini Nvision 40 CrossBeam. 

Atomic force microscope 

Fiber drawing tower Explained schematically in section 13.3.11 For the fabrication of spe- 
cial small hole capillary fibers and customized photonic crystal fibers for metal filling. 

Optics 

Laser Thorlabs HeNe laser HRP 120. Specifications: 12 mW, A = 632.8 nm, linear polariza- 
tion > 500 : 1, mode TEMoo > 99%, 1/e 2 beam diameter 0.88 mm, beam divergence 
0.92 mrad 
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Polarization maintaining fiber Fibercore limited HB600, bowtie structure, single mode for 
A = 633 nm, effective core diameter 3.2 \xm, numeric aperture na = 0.14 — 0.18 

Photo detector (PD1) Thorlabs DET36A High speed si pin photo detector, range A = 
350 — 1100 nm, peak A = 970 nm, active area 3.6 x 3.6 mm, diode capacitance Cd = 
40 pF, responsitivity 7£(633 nm) = 0.4 A/W, bias voltage and max. output voltage 
U = 10 V. 

Amplified photo detector (PD2) Thorlabs PDA100A switchable gain, transimpendance 
amplified si pin photo detector, reverse biased, range A = 400 — 1100 nm, peak 
A = 970 nm, active area 09.8 mm, responsitivity 7^(633 nm) = 0.35 A/W, adjustable 
gain of OdB, 10 dB, 20 dB, 30 dB, 40 dB, 50 dB, 60 dB, 70 dB with the bandwidth 
decreasing for increasing gain from 1.5 MHz to 2 kHz. 

Measurement instruments 

Lock-in amplifier Stanford Research Systems RS830 DSP digital lockin amplifier, external 
inputs with a/d converters, GPIB controlled. 

Lock-in amplifier module 1 Femto LIA-MV-200 single phase digital lockin module, ana- 
log, magnitude signal output. 

Lock-in amplifier module 2 Femto LIA-MVD-200 dual phase digital lockin module, ana- 
log magnitude, phase and quadrature signal output. 

Digital multimeter Agilent 4110A digital multimeter, GPIB controlled and synchronized 
to the RS830. 

Motorized precision rotation stage Physik Instrumente (PI) M-037.DG servo motorized 
, referenced precision rotation stage with PI mercury servo motor controller, usb 
programmable. 
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